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Abstract 

5^ . The magnetic interactions of the two electrons in hehum-hke ions are studied in detail within the 

framework of Relativistic Schrodinger Theory (RST). The general results are used to compute the 
^ ' ground-state interaction energy of some highly-ionized atoms ranging from germanium (Z = 32) 

up to bismuth [Z = 83). When the magnetic interaction energy is added to its electric counterpart 
resulting from the electrostatic approximation, the present RST predictions reach a similar degree 
of precision (relative to the experimental data) as the other theoretical approaches known in the 
literature. However since the RST magnetism is then treated only in lowest-order approximation, 
further improvements of the RST predictions seem possible. 







I. INTRODUCTION AND SURVEY OF RESULTS 



The overwhelming success of relativistic quantum field theory (especially quantum elec- 
trodynamics) could tempt one to belief that any relativistic form of quantum theory must 
necessarily include the virtual degrees of freedom of matter. Indeed it is well-known that 
the virtual processes show up experimentally in form of radiation corrections, vacuum polar- 
ization, pair creation and annihilation etc. 1]. If such a viewpoint would be true, it would 
logically not be possible to construct a consistent relativistic quantum mechanics for inter- 
acting A^— particle systems whose non relativistic description presents no problem at all P| . 
However it seems that there actually exists no convincing argument why such a relativistic 
quantum mechanics for systems of fixed particle number N could not exist. Quite on the 
contrary, since the virtual processes of relativistic field theory (e.g. the Lamb shift) induce 
only small corrections of the results obtained from the first-quantized quantum mechanics 
theories, one concludes that a consistent relativistic quantum mechanics would correctly de- 
scribe the main relativistic effects occurring in the many-particle systems (e.g. heavy atoms) 
but would leave the small corrections due to the virtual processes to the impact of quantum 
field theory. 

In this context, one could now think that the desired relativistic quantum mechanics 
should be deduced as a certain kind of limit process from quantum field theory itself, namely 
through fixing the particle number by neglection of just the virtual processes. Such a deduc- 
tion has been tned already long ago by Bethe and Salpeter BQ but the resnltin. iV-particle 
wave equations are plagued by many deficiencies and interpretation problems p, |6|| . There- 
fore a new approach to relativistic quantum mechanics for many-particle systems has recently 
been established in form of Relativistic Schrodinger Theory (RST) 3, 0, ^, Q| . This ap- 
proach is a relativistic gauge theory based upon the gauge group U(N). The gauge group 
is broken down to the abelian subgroup U(l)xU(l). . . xU(l) which then describes the elec- 
tromagnetic interaction of N particles, whereas the frozen gauge degrees of freedom are 
responsible for the exchange interactions. The non-relativistic limit of this theory coincides 
with the well-known Hartree-Fock approach 

The present paper presents a test of the practical usefulness of RST, namely by elabo- 
rating its predictions for the ground-state interaction energy of the two electrons in highly- 
ionized helium-like atoms where relativistic effects play a non-negligible role. In contrast to 
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the two-particle ground-state energy itself, the corresponding interaction energy is directly 
accessible to the experiments which thus provide sufficient data for testing various the- 
oretical a ppr oaches: relativistic 1/Z expansion multiconfiguration, Dirac-Fock method 



(MCDF) 
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i, 3' relativistic many-body perturbation theory (MBPT) all-order 
technique for relativistic MBPT j]^. These four theoretical approaches have been chosen 
to be opposed to the experimental data; and the present paper adds now the corresponding 
RST results which are found to describe the experimental data with a similar degree of 
precision as these other theoretical approaches (table II). This success is attained by going 
beyond the electrostatic approximation ([20]) and taking into account also the magnetic 
interactions which implies an improvement of the RST predictions by (roughly) one order 
of magnitude. However one should observe here that these RST results have been obtained 
in a preliminary way by a very rough approximation technique, namely by neglecting the 
non-abelian character of RST through linearizing the gauge field equations (Sect. 4B) and 
furthermore by treating the magnetic interactions only in their non- relativistic limit. Thus 
the potentiality of RST is not yet fully exploited and further improvements of the RST 
predictions seem still possible. 

The procedure is as follows: In Sect. II we present a brief survey of RST in order that 
the main arguments of the subsequent discussions can be understood without looking up 
the whole development of the theory in the preceding papers. Here the emphasis is laid 
upon the construction of an energy functional Et, namely as the integral (over all 3-space) 
of an energy density Too(r), see equation ()2.44j) below. Clearly the reason is that, with such 
a functional at hand, the energy level system of the bound A^— particle systems may be 
determined as the value of that functional i?T upon the stationary bound solutions of the 
RST field equations. 

This kind of solutions is then described in great detail in Sect. Ill with an explicit pre- 
sentation of the mass-eigenvalue equations (see equations (I3.4a|l - (j3.4bj) below). Since this 
eigenvalue problem is too difficult to be solved exactly, one has to resort to appropriate 
approximation techniques (as is mostly the case in quantum field theory). Here the logi- 
cal structure of RST suggests to ffist neglect the magnetic interactions between the parti- 
cles ( "electrostatic approximation") which leaves us with a simpler eigenvalue problem (see 
equations ()3.6a|) - ()3.6b|) below). This truncated system is then solved numerically for the 
ground state of the two-particle system in the Coulomb field ^'^^^AQ{r) ()3.9j) where the nu- 
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clear charge numbers (zcx) range from Zex = 32 (germanium) up to z^x = 83 (bismuth). The 
corresponding ground-state interaction energy Ai^^grp (in the electrostatic approximation) 
is then compared to the corresponding experimental values Ai^exp (see table I); and it is 
found that there is a discrepancy between the theoretical (RST) and experimental values 
extending from 1.7 eV for germanium up to 11.5 eV for bismuth. However some intuitive 
arguments indicate that the observed discrepancy should actually be due to the neglection of 
the magnetic forces. When the latter forces are taken into account, there arises a structure 
equation which specifies the desired interaction energy AE in terms of the electromagnetic 
coupling constant (-ZexCts) and two functions and which are only weakly depending 
upon that coupling constant, see equation (|3.59|) and the last two columns of table I. 

In Sect. IV the hypothesis of the magnetic origin of the discrepancy between the elec- 
trostatic RST predictions Ai^^grp and the experimental data Ai^exp is inspected very thor- 
oughly by working out the RST field theory of atomic magnetism. Here it seems reasonable 
to consider first the magnetic effects in the lowest-order approximation, i.e. we linearize the 
non-abelian gauge field equations and additionally we resort to the non-relativistic limit of 
the electrostatic wave functions for calculating the magnetic energy contributions. It is very 
striking to observe that the magnetic and electric contributions display certain dissimilar- 
ities which can however be revealed as necessary consequences of the principle of minimal 
coupling and Lorentz invariance of the theory. Though one obtains a very plausible result 
for the magnetic energy contribution AE^^\ see equation ()4.52|1 below, one nevertheless 
may wish to become convinced by an independent argument which supports our claim that 
magnetism is correctly incorporated RST by the present approach. 

And indeed, such an additional argument in favor of the RST picture of atomic magnetism 
can be supplied, namely by considering the interaction of the bound particles with an external 
magnetic field Hf,^ (Sect. V). Here it can be shown that the interaction energy between the 
particles and the external magnetic field Hex exactly agrees with the conventional results 
for the Zeeman effects; i.e. the magnetic RST energy coincides with the expectation value 
of the conventional Zeeman Hamiltonian, see equation (j5.2j) below. 

Being thus convinced of the physical correctness of the RST picture of atomic magnetism, 
one can return to the hypothesis of the magnetic origin of the numerical gap between the 
electrostatic RST predictions A£^^grp and the experimental data AE^xp of table I. In Sect. VI 
we apply the general results of Sect. IV in order to once more calculate the ground-state 
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interaction energy for the helium-like ions from 2;ex = 32 (germanium) up to z^^ = 83 
(bismuth), but now with inclusion of the magnetic interactions. Here we are satisfied for the 
moment with their linearized description in the non-relativistic limit. Amazingly enough, 
this simple treatment of atomic magnetism is sufficient in order to close the gap between the 
RST predictions AE'j^grp and the experimental data AE'exp up to less than 0.5%, see table II 
in comparison to table I. This must be considered as a rather convincing argument in favor of 
the "magnetic" hypothesis for the observed discrepancy between Ai^^grp and Aii^exp- It is also 
very satisfying to observe that, with the inclusion of the magnetic effects, the corresponding 
RST predictions AE^^^ are now as close to the experimental values Ai^exp as is the case 
with the other approximation methods known in the literature: Relativistic Many-Body 



Perturbation Theory 
Fock method (MCDF) 



17 



14 



15 



1-order technique for MBPT fl^, Multi-Configuration Dirac- 
[3| and relativistic 1/Z expansion [l^, see table II. (The 
predictions of these theoretical approaches have been listed in table III of ref. 1^. Clearly 
one expects that the RST predictions will be further improved by taking into account also 
the relativistic effects of atomic magnetism and retaining the non-linear terms due to the 
non-abelian character of the gauge group U(N) (separate paper). 

Finally it should be stressed that the numerical success of the magnetic hypothesis is 
mainly due to the application of the non-abelian U(2). The reason is that after the "abelian 
symmetry-breaking" (Sect. 2) the frozen gauge degrees of freedom of U(2) imply the existence 
of a "magnetic" exchange vector potential B[r) which plays an analogous part with respect 
to the X- and y-axis as the magnetostatic potentials Aa{r) with respect to the z-axis of the 
coordinate system. Indeed when the spins of the two ground-state electrons are oriented 
along the z-axis, the corresponding spin-spin interaction energy Ae!^'^ (see equation ()6.39|) 
below) is well suited in order to contribute to closing the gap between Aii^^grp and Ai^exp, 
however only one third of the gap could be closed in this way. The other two thirds of the 
gap must be filled by the spin-spin interaction energies Ae!^^ (jOT^ and Ae!^^ (|6lTh|l due 
to the spin orientation along the x- and y-axis (see fig. 1). But the latter interactions are 
described by just that exchange vector potential B{f) being due to the frozen gauge degrees 
of freedom of U(2). Thus each axis of spin orientation contributes the same magnetic energy, 
and this is nothing else than a demonstration of the isotropic geometry of the ground-state. 
In this sense the choice of the non-abelian U(2) is seen to be necessary just in order to 
guarantee the isotropy of the ground state. 
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II. RELATIVISTIC SCHRODINGER THEORY 



In order to let the subsequent elaboration of arguments appear sufficiently self-contained, 
we ffist present a brief sketch of the Relativistic Schrodinger Theory (RST), which places 
emphasis on a closer inspection of the energy functional, since this will ultimately yield the 
energy level system of the bound many -particle systems. For more details and deductions, 
the interested reader is referred to the preceding papers, e.g. refs. f j7|)-([ll|). 



A. RST Dynamics 

As for any field theory of matter, the fundamentals of RST consist in a basic system 
of field equations, which is subdivided into three coupled subystems: matter dynamics, 
Hamiltonian dynamics and gauge field dynamics. 

(i) Matter Dynamics 

When matter occurs in form of pure state its distribution over space -time is governed 
by the Relativistic Schrodinger Equation ( RSE ) 

ihcv^^i = n^^ . (2.1) 

If it is more adequate to describe matter by a mixture, being characterized by an intensity 
matrix T, one requires I to obey the Relativistic von Neumann Equation ( RNE ) 

v^i = ^^[i-n^-n^-i] . (2.2) 

Clearly, a pure state \1/ can be considered as a special type of mixture, namely that for which 
the intensity matrix degenerates to the tensor product of \E' and its Hermitian conjugate ^, 
i.e. 

J^^®i|r. (2.3) 

In the present paper we are restricting ourselves to the investigation of stationary bound 
systems being described by pure states. Since a pure state \1/ of an A^-fermion system has 
4:N components, we have to consider a C^^ - realization of RST. 

(ii) Hamiltonian Dynamics 

The Hamiltonian occuring in the matter field equations f 12.11 ) and f l2.2l ). is itself a 
dynamical object obeying its own field equations, namely the integrability condition 

V^H, - V,n^ + ^ [n^,, Hu] = ihcJ^^, (2.4) 
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and the conservation equation which reads for fermions 



V^'H^ - -n^ .n^ = ^hc^—j ■ 1 - thcJl^^J^^, . (2.5) 

For bosons the spin term (last term on the right-hand side) is omitted. The meaning of 
the integrabihty condition f l2.4l ) is to ensure the vahdity of the bundle identities 

[V^V, - V,V^] ^ = J-^, . ^ (2.6a) 
[D^V, - V,V^] 1 = \T^,, T] . (2.6b) 

Similarly, the conservation equation ( 12.51 ) guarantees the existence of certain conservation 
laws, e.g. the charge conservation 

V% = , (2.7) 
or the energy - momentum conservation for free particles (Dirac particles here) 

V'^(^>r^. = , (2.8) 

whose energy- momentum density is denoted by ^^^f_iu- 

Clearly, if the considered A^- particle system is interacting with an external field ^^^^J-'^u 

(ex)T7 _ _ • (ex)p . 1 — f _ (S)-r- (n q^ 

the energy - momentum density ^^^^u cannot obey the conservation law ( 12.81 ) without being 
complemented by the energy- momentum density ^^^^^u due to the external interaction. 
Consequently, the conservation law ( 12.80 has to be replaced by the balance equation 

y,. ((S>^^^ ^ (es>^^^) ^ _(xe)^^ ^ (2.10) 

where the Lorentz - force density ^^'^^fu is given by 

(^^)/, = -hcF^^^'^Y ■ (2.11) 

Here, ^'^^^j^ is the external four -current generating the external field ^'^^^F^^, according to 
Maxwell's equations 



VMex)^^^ = 47ras(''"^J. (2.12) 



as 
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and F^y is the coherent electromagnetic field strength generated by all the particles of the 
system: 

F,^ = J^^^{^'^^,A ■ (2-13) 

This means that the source of this total field F^i, is the total current of the A^- particle 
system 

V^F^, = -Anasj, , (2.14) 

which obeys the strict conservation law f l2.7l ).If one prefers to think of a closed system, one 
adds the energy - momentum density '^^^^^i, of the external source 

yM(ex)j.^^ = ("^)/, (2.15) 

and then finds the desired closedness relation by combining equations f 12. 101 ) and f 12. 151 ) 

((S)r^^ + (cs)r^^ + (cx)j.^^) = . (2.16) 

Finally, let us also remark that the conservation equation ( 12.51 ) is equivalent to the 
following condition upon the Hamiltonian: 

V^" -71^ = Mc^l . (2.17) 

Here, the total velocity operator generates a (4A^) - dimensional representation of the 
Clifford algebra C(l, 3), i.e. 

F^-F, + F, ■F^ = 2(?^,-l, (2.18) 

with the corresponding generators H^^y of the group Spin{l, 3) being given by 

^,. = \[T„T,] , (2.19) 

see the right-hand side of the conservation equation 12.51 The physical meaning of F^ is to 
build up the total current according to 

j^ = tr{J-F4 , (2.20) 

i.e. for a pure state f l2.3l ) 

= ■ F^ ■ ^ . (2.21) 



Since by means of the latter form f l2.17l ) of the conservation equation the RSE f 12.11 ) can 
be converted to the A^- particle Dirac equation 



ihT^'Vf,'^ = Mc^ , (2.22) 

it thus becomes a simple matter to verify directly the charge conservation law f l2.7l ). This 
form of the charge conservation does not only work for pure states but equally well for mix- 
tures, which can be proven by resorting to the RNE ( 12.21 ) instead of the RSE ( 12.11 ). Observe 
also that the non - relativistic limit of the iV- particle Dirac equation ( 12.221 ) coincides with 
the well-known Hartree-Fock equations jll| . 
(iii) Gauge Field Dynamics 

In order to have a closed dynamical system, one finally has to specify a field equation for 
the bundle curvature ( "field strength") J^^^. This is the non-abelian Maxwell equation 

V'T^, = -AmasJu • (2.23) 

Similarly, as it was done for the field strenght JF^^/ ( 12.91 ). the current operator J'^ may also 
be split into an external and a system part 

= ^"^^M + ^'^^M (2-24a) 
(ex)^^ = ■ 1 (2.24b) 

so that the non-abelian Maxwell equation f 12. 231 ) can be required to decay analogously into 
two parts 

p/^Mj-^, = -47rz«g(^^)j, (2.25a) 
V''^'^^, = -47rz«s^^^X . (2.25b) 

Here the external part f l2.25an gives the former ( abelian ) Maxwell equation f 12. 121 ). whereas 
the "total" Maxwell equation f 12. 141 ) emerges as the trace part of f l2.25bl ) with 

i/. = -^^tr{J-4 . (2.26) 

The essential point with the gauge structure of RST refers to the right choice of the gauge 
group. In conventional quantum ( field ) theory, the group U{1) is evoked in order to describe 
the electromagnetic interactions. However, in RST the - particle bundle of wave functions 



8 



\I/ is the Whitney sum of the corresponding single - particle bundles which suggests to adopt 
the A^-fold product group f/(l) x f/(l) x . . . x f/(l) as the adequate subgroup of U{N) for 
the description of the RST gauge interactions. Though the electromagnetic interactions are 
thus adequately incorporated into RST by means of this "ahelian symmetry breaking", one 
needs the residual structure of the embedding group U (N) in order to take into account also 
the exchange interactions between the particles. 

Considering for instance the two-fermion case ( = 2), one has to deal with two "elec- 
tromagnetic generators" (a = 1,2) and two "exchange generators" X:X iii order to span 
the gauge algebra u(2). The bundle connection ( "gauge potential" ) may then again be 
splitted up into an external (ex) and internal system part (S) 

A, = (^^'U^ + (^U^ (2.27a) 
(^^U^ = -^(^^U^ . 1 , (2.27b) 

where the internal part '•'^M^ decomposes with respect to the u(2) basis {tq, x, x} as follows 

^'U, = A\Ta + B,x-B;x. (2.28) 
Clearly, the connection Afj, generates its curvature jF^j^ as usual 

~ ^ iiA^ -\- \Afi., Au] , (2.29) 

and by decomposing the internal part ^^"'T^y ( 12.91 ) in a similar way 

= F^^^T, + G,yx - G%x , (2.30) 

the curvature components -F^^j,, G^y are found in terms of the connection components 
A^^, as follows: 

^'""K. = V^(^^U, - V,^'^^ (2.31a) 

F\y = V,A\ - VyA\ + I [B,B: - B^B;] (2.31b) 

= V,A\ - VyA\ - I [B,B: - ByB;] (2.31c) 

= V,By - VyB, + I [A\ - A\] -By -I [A\ - A\] ■ B, (2.31d) 

G% = V,B: - V.B; - I [A\ - A\] -Bl + i [A\ - A\] • B; . (2.31e) 
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The total field strength F^i, f 12. 131 ) appears now as the sum of the "electromagnetic" field 
strengths F""^^ (ESIEl) - (|2ISIcl) , i.e. 

F,, = F\^ + F% , (2.32) 

and therefore F^^j is generated by the "total potential" 

A, = A\ + A\ (2.33) 

in a strictly abelian way 

F^, = V^A, - V,A^ . (2.34) 

Recalling also the abelian character of the "total" Maxwell equations f 12. 141 ) we see that the 
total objects of the system A^,F^u,jfj, obey a strictly abelian structure due to the "total" 
group U{1), just as it is the case in conventional (quantum) electrodynamics. 

In contrast to these "total" properties, describing the - particle system as a whole, the 
"internal" gauge structure is truly non- abelian on account of the exchange interactions. 
More concretely, when the current operator '•^^JT'^ f l2.24al ) is decomposed as 

= ^ {-k\Ta + hlx - Kx] , (2.35) 

the "internal" Maxwell equations f l2.25bl ) read in component form 

WF\^ + t [B^Gl, - = -A7^ask\ (2.36a) 

V^F^ - I [B^G% - B*^G,,] = -Anask\ (2.36b) 

WG^, + I [A^^ - A^^] G^, - iB^ [F\^ - F%] = A-Ka^K (2.36c) 

V^G;, - z [A'^ - A^^\ G% + iB*^ {F\, - F%] = Ai^a^K ■ (2.36d) 

Adding here both equations ( I2.36al ) and ( I2.36bl ) yields just the "total" Maxwell equations 
( 12.141 ) with the total current now being found in terms of the electromagnetic single - 
particle currents k""^ as 

J, = k\ + k\ . (2.37) 

Finally, one may also consider the non- abelian "charge conservation". Observe here, 
that the curvature J-'^^ must obey the following bundle identity (valid in any space -time 
with vanishing torsion ) 

VVJ^^^ = . (2.38) 
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Therefore differentiate once more the internal Maxwell equations f l2.25bn and find the fol- 
lowing source equation in operator form 

V^'^^^ = . (2.39) 

Written in component form, this source relation reads 

V^A;^^ = i [B^'h^ - B*''h*^] (2.40a) 

V^'k\ = -i [B^h^ - (2.40b) 

V^hf, = I {A^^ - A^^\ + iS*^ [fci^ - k^^ (2.40c) 

V^/i; = -I {A^^ - A^^\ h; - iB'^ [k\ - k\] . (2.40d) 

As a consistency test, one may add up the first two source relations f l2.40an - f l2.40bn in 
order to find just the former charge conservation law ( 12.71 ) for the total current f 12. 3 71 ). 

Clearly, the crucial point with this non-abelian gauge structure of RST is now whether 
one can identify certain experimental facts which support the non-abelian hypothesis. Since 
the latter aims at the exchange interactions of many -particle systems, one may consider the 
stationary states of bound A^- particle systems (e.g. many- electron atoms) which occur in 
form of entangled states implying the emergence of "exchange energy" . Or in other words, 
verification of the RST predictions may be obtained by inspection of the atomic energy level 
systems, which are determined to a considerable extent by the exchange interactions. In 
order to prepare such a test of RST, it is very instructive to consider first an appropriate 
energy functional (-Et, say) for the RST field configurations. 



B. Energy Functional 

Within the context of a relativistic field theory, one would like to define the energy E^y 
as the space integral over the time component '•^^oo of an appropriate energy - momentum 
tensor ^^'^^u- Equation ( 12.161 ) suggests to take for ^^^^u the sum of the internal part ^^^^v 
of the system and the interaction part ^^^^^u with respect to the external source 

^^Ku = ^'K^ + ^-Ku ■ (2.41) 

Since every A^- particle system is composed of matter and gauge fields, mediating the inter- 
actions of matter, one will build up the internal density ^^XT^^ by a Dirac matter part ^^^^u 
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and a gauge part ^^XT^^: 

(^>r,. = (°>r,. + (^>r,. . (2.42) 

Thus the total density ^^^fj,u f 12.411 ) of the - particle system consists of three parts 

^^K^ = + («>r,. + (->r,. . (2.43) 

Consequently, the total energy 

Et = j d¥(T>roo(r-) (2.44) 
will also be built up by three constituents 

E^ = Eo + Eg + Ees , (2.45) 

with the self-evident definitions 

E^ = j d¥(°>roo(f) (2.46a) 

Eg = J d¥(«>roo(f) (2.46b) 

Ees = [ dV-'(^^>roo(r) . (2.46c) 



Here the simplest part is E^g f l2.46cn . because the corresponding energy - momentum 
density ^'^^^^u is in a bilinear way composed of the the external field strength ^^^^F^i, ( 12.91 ) 
and the total field strength F^^ f 12. 131 ): 

^"'X. = l^'^^^i^^.A • + F^, • - i(7^.(^^)F.A ■ F-^^l . (2.47) 

Consequently, the external interaction energy E^s ( I2.46cl ) is given by 

he 



E,, 



Aira 



J d¥ l^e. ■ F + i?ex ■ i?} = Eit^ + Ei^^ , (2.48) 



where the three- vectors E{r) = {E^} and H{r) = {H^ of electric and magnetic field 
strengths are given in terms of the components of the field strength tensor F^^y ( 12.321 ) by 

E^ = Foj (2.49a) 
if^ = V,F,',etc. (2.49b) 

— * — * 

If the external fields E^^, H^x are adopted to be constant (as it is usually done for treating 
the Zeeman and Stark effects), the external energy E^s ( 12.481 ) becomes infinite and one has 
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to substract the infinite contribution due to the external source generating the homogeneous 
fields Ec^, Hex (see the discussion below equations f 13. 341 ) and f l5.8l ) ). 

Next, considering the gauge field density one finds this object to appear as the 

difference of the electromagnetic energy - momentum density ^^'T^^ and the exchange density 

(^>r,. = (^>r,. - (c>r,. (2.50) 

with the energy -momentum content of the electromagnetic modes being given by 

'""^u = {f^FV + F\,FV - \9,.F\,F-^>^'^ , (2.51) 

and similarly for the exchange modes of the gauge field system 

^'^^fj.u = -^^^^ |g'^aG'*'^ + G'*^G/ - -g^ivGlxG''^^ ■ (2.52) 

Naturally, according to this subdivision, the gauge field energy Eq f l2.46bn is also split up 
into two parts 

EG = En-Ec, (2.53) 
where the electromagnetic energy is given by 

E^ = I dV-'(^>roo(f) = ^ 1 dV- [e, ■E, + H,- #2} , (2.54) 

and similarly for the exchange energy Eq 

Ec = j d¥(c>roo(r) = J d¥ |x* ■ X + F* • f} . (2.55) 

Clearly, the "electric" and "magnetic" exchange three - vectors X = {X^} and Y = {Y^} 
are defined analogously to their electromagnetic counterparts f l2.49an - f l2.49bn 

X^ = Goj (2.56a) 
= y^'Gu . (2.56b) 

Now it is evident that both the electromagnetic energy E^i ( 12.541 ) and the exchange 
energy Eq f 12. 551 ) may be split up further into their "electric" and "magnetic" constituents, 
i.e. we put 

En = 4'^ + 4"^ (2-57a) 
Ec = 4'^ + 4^ , (2-57b) 
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with the "electric" parts being defined in an obvious way through 

n(e) ^ /■ 3 



47ras 

and similarly for the "magnetic" parts 



= / dVEi ■ E2 (2.58a) 

47ras J 

= / d'fX* ■ X , (2.58b) 



Anas J 



E)~f> = / dVr* ■ Y . (2.59b) 

The important point with this subdivision into "electric" and "magnetic" parts of the gauge 
energy is, that the "electric" contributions E^j^ and -E^^ will turn out to be much larger 
than their "magnetic" counterparts E^^^ and E^^\ so that the latter ones may be treated 
as small perturbations. The main effect is therefore due to the "electric" contributions 
and represents thus by itself an acceptable first approximation to the experimental data 
{Electrostatic Approximation). In this sense we may rearrange the gauge field energy Eq 
(E3S1) as 

= 4"^ - 4'^ + (4"^^ - 4^) , (2.60) 



so that the "magnetic" contributions (in brackets) can be omitted for the electrostatic 
approximation. This approximation is then based solely on the truncated gauge field energy 



£;(eh) 

G 



Egi^) = - E^ (2.61) 



and thus misses the magnetic energy E^™^^ 

= - E^f) , (2.62) 

which will be considered as the "magnetic correction" of the electrostatic energy 

Finally, the matter energy Ed f l2.46al ) has to be considered. The corresponding energy - 
momentum density '^^^^j.u carried by the wave function \Ef is given by 

^''%u = ^ {^^r^ {V,m) - {V,^>) V^m + ^>T, {V^^) - (P^fr) r,^} . (2.63) 

Intuitively, one expects the associated matter energy Ed to consist of the rest mass energy 
2Mc^ of both particles ( take for the moment N = 2) plus their kinetic energy E^in, i-e. one 
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expects 

2 

Ej,=2Mc^ + J2EMn(a) . (2.64) 



a=l 



Indeed, this result can actually be deduced in the non-relativistic limit j21[. In order to 
have a brief sketch of the proof, eliminate the time derivative X'o^ from the matter energy 
density ^^^oo 

(°>roo = ^ {^ro (Vo^) - {Vo^) To^} (2.65) 
by means of the Dirac equation ( 12.221 ) and find 

(°>roo = Mc^^^ + ^ { (Vj^) r% - ^PPj^} . (2.66) 

Integrating now over all three -space and carrying out the non-relativistic limit in a con- 
sistent way yields just the claimed result ( 12.641 ). where the correspondence of rest mass 
terms and kinetic energy terms should be obvious by comparing both equations f 12. 641 ) and 

(HH). 

Observe also, that for the stationary bound field configurations ( as solutions of the mass 
eigenvalue equations) the mass eigenvalues should appear somewhere in the energy func- 
tional. For one -particle systems, one may even expect that the field energy Eq^ has to 
coincide with the mass -energy eigenvalue (M*c^, say) of the bound particle [lfll |: 



Ft = M,c^ . (2.67) 

The relationship between the field energy i^x and the mass eigenvalues for few -particle sys- 
tems has now to be studied in more detail, together with their dependence on the "electric" 
and "magnetic" interactions. 



III. ELECTROSTATIC APPROXIMATION 



The relative magnitude of the "electric" and "magnetic" interparticle interactions in a 
stationary bound system strongly depends on the value of the electromagnetic coupling 
constant as- In view of its smallness (as ~ liy)' "electric" interactions clearly 
dominate the "magnetic" interactions, because their relative magnitude is of the order 
(zex«s)^; where Zex is the number of (positive) charge units located at the center of the 
binding force, e.g. the nucleus j^]. Therefore, at least for the light atoms {zex < 10, ...,20, 
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say), one may speculate that the complete neglection of the interactions of the "magnetic" 
type ("electrostatic approximation" ) yields a first useful estimate of the relativistic energy 
levels of the bound system. The reason for this expectation is that, as a consequence of 
their dominance over the magnetic interactions, the electric interactions will also contribute 
the main part of the relativistic effects. If this presumption is true (to be verified readily), 
one will not only be able to treat the magnetic interactions as a small perturbation but 
furthermore the non- relativistic limit of this perturbation will be all that is needed for 
taking into account the magnetic effect to first order. In chapter VI we will demonstrate 
that this apparently rough approximation scheme gets the RST predictions closer to the ex- 



perimental data of ref 



12j ending up with less than 0, 5% of relative deviation ( see table II ). 



Thus we will proceed as follows: 

(i) elaborating in this chapter the relativistic stationary field configurations as far as they 
are needed for the electrostatic approximation 

and 

(ii) computing the magnetic corrections in the non - relativistic limit in the subsequent 
chapters. 



A. Stationary States 



The stationary bound states are not completely time- independent but have such a tem- 
poral behaviour that the observable quantities of the theory ( i.e. eigenvalues and densities 
of charge, current, energy etc. ) become truly time- independent. Thus the wave functions 
of the considered two - particle system ( = 2 ) vary with time as 



i!i{r,t) 
^2(r,t) 



exp 



exp 



— z- 



h 

f2C 



Moc^ 



h 



-t 



(3.1a) 
(3.1b) 



where Ma (a = 1,2) denote the mass eigenvalues. This stationary form of the wave func- 
tions implies the time - independence of the electromagnetic objects, i.e. the currents /c°^. 
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potentials and field strengths Ea, Ha- 

{r4^{(%(r);4(r) 
{A»J^{('^)Ao(f);L(f)} 
Ea Ea{f^ 

Ha Ha{^ , 

whereas their exchange counterparts adopt the following time dependence: 



(3.2a) 
(3.2b) 
(3.2c) 
(3.2d) 



h^{r,t) = exp 
h*^ir,t) = exp 
B^{f,t) = exp 
B*^{r,t) = exp 
X(r, t) = exp 
Y{f, t) = exp 



'.Ml - Ms 2/ 

I c t 

n 



M1-M2 2, 

—I c t 

n 

.M1-M2 2, 

—I c t 

n 

■.Ml -M2 2, 

2 C t 

n 



Ml -M2 2, 

—2 C t 

n 

.M1-M2 2, 

—z c t 

n 



■ K{f^ 

■ xi^ 

■Y{r) . 



(3.3a) 
(3.3b) 
(3.3c) 
(3.3d) 
(3.3e) 
(3.3f) 



The relativistic eigenvalue problem for the determination of the mass eigenvalues Ma is 
obtained now by substituting the stationary form of the wave functions into the two - particle 
Dirac equation ( 12.221 ) which yields the following coupled system for the spatial parts ipai^ 
of the stationary wave functions ipa{r,t) ( I3.1al ) - ( I3.1bl ): 

i^m.M^ + P'^"^^o(r) + (')Ao(r-)] 70^1 (f) + 5o(f)7o^2(r) = ^ ~ c^i(f) (3.4a) 
^7^ro,^2(r-) + [(^^Uo(r-) + (^Uo(f)] joM^ + B*{f)joM^ = ^ ' . (3.4b) 

Here the time derivatives of the wave functions have brought in the mass eigenvalues Ma, 
and the remaining spatial derivations Dj are defined as follows: 



Mr) = djMr) - i P'%(r0 + ^%(r)] Mr) - iSj(f)^2(r) 



(3.5a) 
(3.5b) 



Now, by its very definition, the electrostatic approximation consists in neglecting the inter- 
actions of the "magnetic" type, i.e. one puts to zero all three -vector potentials Aa{r), B{f) 
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and thereby ends up with the following truncated system: 

^7^5,V^i(f) + [(^^Uo(f) + (2Uo(f)] 70^1 (f) + 5o(r)7o^2(r) = ^ ' cV^i(f) (3.6a) 
^7^a,V'2(r) + [(^^Uo(r) + ('Uo(f)] 7o^2(r-) + 5*(r)7o^i(r) = ^ c V'2(r) . (3.6b) 

Logically, both Dirac spinors ■i/'a (r) ( a = 1 , 2 ) do couple here exclusively to the remaining 
time components '•"^74o(r'), Bo^f) of the corresponding four- vector potentials A"-^, B^. There- 
fore we have to complement the RST - Dirac system f l3.6an - ( I3.6bn merely by the Poisson 
equations for these time components, which of course have to be deduced from the Maxwell 
equations f l2.36al ) - f l2.36dn by neglecting again the spatial components of the four- vector 
potentials: 

A('^)/lo(f) = 47ras^%(r) (3.7a) 
A5o(f) = -A'KasKir) . (3.7b) 

Observing here that the currents fc"^, /i^ are generated by the wave functions through 

k\{f^ = ^a(r)7^V'a(r) (3.8a) 
/i^(r-)=V^i(f)7^^2(r), (3.8b) 

one arrives at a coupled but closed Dirac - Poisson system f l3.6an - f l3.7bn whose solutions 
'?/'a(r'), Ma| constitute what we consider to be the "electrostatic approximation". 
Subsequently, we want to carry out a test of the usefulness of RST by opposing its 
theoretical predictions for the two -particle ground -state in the Coulomb potential 

(^^Uo(f) = Ze.^ (3.9) 
\r\ 

to the experimental data. Since the ground- state is a member of the para-sytem with the 
highest possible symmetry, we try the following ansatz for the Dirac spinors ipai"^)'- 

( ("U+(f)\ 

with the first particle having spin up, i.e. the two component Pauli spinors are chosen as 



1 1 



«0+(r) = (')i?+(r)-CV (3.11a) 



1 1 



(lU (f) = -i (i)i?_(r) ■ C? ■ (3.11b) 
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Here the one - particle eigenspinors of the angular momentum operators obey the rela- 
tions Q, Q 







(3.12a) 


CT 




(3.12b) 




= i{i + i)ff ■ CT 


(3.12c) 


^2 ^j,m 




(3.12d) 



with the composition law for spin - ^ particles 

3=l±\. (3.13) 

Since the two spins are anti - parallel in the ground- state, we put for the second particle 

(2)0^ (f) = (2)i?^(r) ■ (3.14a) 
(2)0_(f) = -i (2)i?_(r) ■ . (3.14b) 

Now, for the ground state, it is reasonable to assume that the radial parts ^'^'>R±{r) of 
both wave functions ipai'f^ are identical 

^^^R+{r) = ^"^^R+ir) = R+{r) (3.15a) 
(r) = (r) = R-{r), (3. 15b) 

and similarly for the mass eigenvalues 

Mi = M2 = Mo . (3.16) 

Thus the original RST eigenvalue system f l3.6al ) - f l3.6bl ) becomes simplified to two coupled 
equations for the two radial ansatz functions R±{r): 

+ [(-Uo(r) + Mr)] . R.ir) = , . R_^r) (3.17a) 

+ lR_^r) - [(-Uo(r) + Ao{r)] ■ i?+(r) = ^° ~ ^ c ■ i?+(r) . (3.17b) 
ar r n 

Observe here that for the para - ansatz f l3.14al ) - f l3.14bn the time component hQ{r) of the 
exchange current f l3.8bn vanishes ( /io(^) = ) so that one can take a vanishing exchange 
potential (-Bo(^) = 0) as solution of the Poisson equation ( I3.7bl ). Therefore the radial 

19 



functions R±{r) do couple only to the electrostatic potentials ^"-^AQ^r) which of course now 
have to be identical, too 

(%(r) ^ (2)^o(r) = Ao(r) (3.18) 
and obey the Poisson equation f l3.7al ) 

AAoir) = Anaskoir) , (3.19) 

with the coinciding electric charge densities ^""^kQ^r) f l3.2al ) being given by 

n2 , /r, / \\2 



(1), / A (2), / A . , / N (R+ir)) + {R-{r)y 



(3.20) 



Thus the formal solution of the Poisson equation ( 13.191 ) is found to be 

Mr) = -^fi^-i^±qL±(^, (3.21) 

47r J \r — r'\ 

B. Non - relativistic Limit 

Since the electrostatic approximation is required to take fully into account the relativistic 
effects, the coupled system ( I3.17al ) - ( 13.201 ) has to be solved numerically and then yields the 
basis for the RST results AE^^^j, (third row of table I). However, for the calculation of the 
magnetic corrections we may restrict ourselves to the lowest -order approximation of the 
solutions ilJaif)- This means that we 

(i) neglect the interaction between both particles (i.e. putting Ao(r) ^ 0) 
and 

(ii) we are satisfied with the non - relativistic approximation of the corresponding solutions 
R±{r). 

Now, by virtue of the first assumption, the interactive system f l3.17al ) - f l3.17bl ) becomes 
decoupled and truncated to the simpler form 

^ + (-)Ao(r) . R.{r) = - ^- + ^ c ■ R.{r) (3.22a) 
dr h 

+ -R-{r) - (-Mo(r) . R^{r) = c ■ R^{r) , (3.22b) 

dr r h 
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whose solutions are given by 



2i 



25| 



R_(r) = A^.v^M-M.r^'exp 



Ob 
Ob 



(3.23a) 
(3.23b) 



Ob 



Me2 



. . . Bohr radius ) . 



Here the parameter u is given by 



ly = -1 + \ I - (zexas 



(3.24) 



Furthermore the mass eigenvalue Mq of the coupled two - particle system f l3.17an - f l3.17bl ) 
degenerates to the well-known one -particle result M* 24 l25l| 



M,=M\ l-{zexas] 



(3.25) 



and finally the normalization constant A^* is computed via the relativistic normalization 
condition 

(3.26) 



j d^fko{r) = 1 



as 



1 V OB 



3+2u 



2MT{3 + 2u) • ^^-^^^ 
The second assumption (of non- relativistic approximation) says that the radial function 

i?_(r) f l3.23bn is neglected against R+{r) and for the latter function ( I3.23al ) one takes its 

o 

non - relativistic limit i?+(r), i.e. 



i?+(r)^i?+(r)=2i 



exp 



Ob 



(3.28) 



which obeys the non -relativistic version of the normalization condition f 13. 261 ) 



dr I R^(r 



1 . 



(3.29) 



Clearly, this non - relativistic limit R+{r) is the ground -state solution of the ordinary 
Schrodinger equation 

^2 



2M 



(3.30) 



2aB 
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and therefore it may appear somewhat amazing that the ordinary Schrodinger equation is 
sufficient in order to calculate the magnetic corrections in the lowest - order approximation. 
However, as we shall readily see, relativity nevertheless enters the "magnetic" exchange 
energy, namely in form of the exchange current h{r) which is built up by the radial functions 
^"•^R±(r) through the combination 

M+(r) = (i)i?+(r) ■ (2)i?_(r) + (2)i?^(r) ■ ^^^R^r) , (3.31) 

see equation f 16. 261 ) below. Observing here the coincidence of the radial functions (")i?±(r) 
for the two -particle ground -state f l3.15an - f l3.15bn . one deduces from the non - relativistic 
approximations ( I3.23al ) - ( I3.23bl ) the corresponding approximation for R+(r) ( 13.311 ) as 

3 r 

(3.32) 



(r) 4 (zextts) -T- exp 



2 ^ex"^ 

Ob 



C. Mass Eigenvalues Ma and Field Energy 

Intuitively, one should suppose that the total energy i?T n2.45l )is related to the mass 
eigenvalues Ma in some way. Indeed, this desired relation can easily be obtained by inserting 
the stationary wave functions ipaif) ( see the right - hand sides of f 13.1 ah - f l3.1bn ), potentials 
A'^^ f l3.2bn and f l3.3cl ) - f l3.3dn into the matter density ^^-Too f 12. 651 ) and integrating over 
all three -space. In this way, one finds quite generally the following result for the matter 
energy Ed ( I2.46al ) in terms of the mass eigenvalues Ma 

Ed = zi- Mic^ + Z2 ■ Msc' - E^^^ (3.33) 
+ hcj d¥ {(^Uo(f) ■ (%(f) + (2)^0 (f) ■ (^)A;o(f)} 

+ hcj d¥ {Eo(f) ■ /io(r) + B*{r} ■ h;{r)} . 
Here it is presumed that the external interaction energy of the electric type E^t^ (E3H1) 

can 

be recasted into the following form by means of GauB' integral theorem and the Poisson 
equations ( I3.7al ): 

Eil^ -he J d¥('^"Uo(f) {^%ir) + (%(f)} . (3.34) 

This presumption implies that the external field Ef,j.{r) is well- localized; if this is not true 
(e.g. by considering a homogeneous field E^x), the GauB' surface term becomes infinite at 
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infinity and has to be omitted. The reason is that such an infinite surface term has to be 
considered as part of the (infinite) energy of the external non-locahzed source generating 
the homogeneous field E^r^. Moreover, the former normalization condition f 13.261 ) upon 
the wave functions ipai^, being adequate for the electrostatic approximation, has to be 
generalized in the presence of magnetic interactions to 

d¥(%(f)=z, , (3.35) 

unity (for the relativistic normalization of wave 



with the real numbers Za still close to 



functions in the general case see ref . j20|, |2l[ ). 

Now one substitutes the present result for the matter energy Ed f 13. 331 ) into the total 
field energy E't f 12. 451 ) which then reappears in the following form: 

Et = zi- M,c^ + Z2 ■ M^c^ + ^ir^ (3-36) 
+ E^^^ + hcj d¥ {(^Uo(f) ■ (')fco(r) + ('Uo(f) ■ ^%{f)} 

- + hcj d¥ {i?o(f) ■ hoir^ + B*ir) ■ /i*(f)} 

+ 4-) - 4) . 

However, the gauge objects emerging here (namely the potentials *^"Mo(r) and the charge 
densities ^"-^kQ^r) ) are not independent, but rather are connected by the Poisson equations 
f l3.7an : and this establishes the following relationship among the electrostatic energy con- 
tributions 



2 _ 

This result again comes about by applying GauB' integral theorem to the internal energy 
functional E^^^ f l2.58al ) and using also the Poisson equations f l3.7al ). Clearly, there exists 
an analogous relationship for the "electric" exchange energy 

4^^ (ElHEl), namely 

^(h) f ^^^^^ . f^^^^ ^ . f^*^^^ ^ f ^3^^^^ . ^.^^ ^ ^338) 

2 J 47rasaM J 

which will be discussed in greater detail below. By means of these relationships the total 

energy Et ( 13.361 ) can now be simplified to the plausible result 

Et = zi- Mic^ + Z2 ■ Mac' + 4"^^ - a4''^ + 4"''^ • (3-39) 
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Here the "magnetic" contribution E^q^^^ f 12. 621 ) reads in terms of the vector potentials Aa{r) 
and i?(r) 

^(mg) ^ _K f ^3- . ^^(^ + A?2(f) . Ai{r) + h{r) ■ B{r) + /I*(f)5*(f)|(3.40) 



2 

he 



and similarly its "electric" counterpart is defined through 

.(ch) ^ j-,{eh) he 



AE^ = Er + ^— / d'rB\r) ■ B{f) (3.41) 
= / d¥ {(lUo(f) ■ (%(r) + (2Uo(f) ■ (i)fco(r) + i?o(r) ■ h,{r) + i^o* (r ) /i* (f ) } 



2 „ 

+ „ /" dV5*(r) ■ 5(r) . 

47rasaM ^ 

This now is a very interesting result, because it provides us with much insight into the 
difi^erence between the "electric" and the "magnetic" interactions. First observe that either 
of the two mass eigenvalues (a = 1,2) already completely incorporates the "electric" 
plus "magnetic" interactions of one particle with the other, respectively. This can be most 
clearly seen from the original form of the eigenvalue equations f l3.4an - f l3.4bn . namely by 
multiplying those equations by ipi or ?/'2, resp., from the left and integrating over all three - 
space and finally resolving for the eigenvalues Ma- Therefore the sum of the mass -energy 
eigenvalues Za ■ MaC^, occurring in the total energy E^ ( 13.391 ). counts the electromagnetic 
interaction energy between both particles twice whereby the "electric" parts enter with a 
positive sign and the "magnetic" parts with a negative sign. Now, one of the pleasant 
features of the RST energy functional E^ f 13. 391 ) is that this double - counting is corrected 
by subtraction of the "electric" interaction part AEq^^ and addition of the "magnetic" 
interaction part i^^^^The relatively good agreement of this theoretical picture with the 
experimental data ( see table I and II ) may be taken as a support of both the RST philosophy 
and of its practical consequences which will be considered readily. 

As far as the method of electrostatic approximation is concerned, it is self - suggesting 
to omit the magnetic terms in the total energy E't ( 13.391 ). quite similarly to the way the 
electrostatic eigenvalue system ( I3.(ial ) - ( I3.6bl ) was obtained from the general sytem ( I3.4al ) - 



( I3.4bl ) by omitting the "magnetic" potentials Aa{r) and B{f). This leads us to define the 
electrostatic approximation Ej- of the total -Et by 

Et = M^c^ + Mac^ - ^E^G^^ , (3.42) 
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where the mass eigenvalues Ma are the eigenvalues of the truncated system f l3.6al ) - f l3.6bl ) 
and the "electric" interaction functional Ai?^^'' is to be taken upon the solutions ■?/'a(^) of 
this truncated system: 

/\Ef^ ^ -y y dV-'{(%(f) • (%(f) + (2)Ao(f) ■ (%(f) + i?o(f) ■ /io(f) + 5°*(f) • hl{r)] . 

(3.43) 

The corresponding numerical results for the ground -state interaction energy Ai^^grp are 
presented in the third column of table I. 

D. Numerical Results 

Within the framework of the electrostatic approximation all prerequisites are now at hand 
in order to test the quality of this type of approximation. The numerical solution of the 
electrostatic eigenvalue system ( I3.17al ) - ( I3.17bl ) provides us with the radial ground -state 
functions R±{r) together with the mass eigenvalue Mq. Consequently, these results may 
now be substituted into the electrostatic energy functional Et ( 13.421 ). whose value upon the 
calculated ground- state solutions |-R±(r);Mo| appears then as 

Et 2MoC^ + nc j d¥Ao(r) • fco(r) (3.44) 

where the electrostatic potential v4o(r) and the charge density fco(^) are given by equations 
f l3.20n - f l3.2in . Observe here also that the exchange potential BQ{r) as a solution of the 
Poisson equation ( I3.7bl ) vanishes because the exchange current hfj_{r) ( I3.8bl ) has a vanishing 
time- component {ho{r) = 0) for the stationary field configurations being defined through 
equations (EHil) - (ElM) and (Ellil) - (EHB). 

Now the problem is how to judge of the present RST prediction E^^ f 13. 441 ) for the two- 
electron ground -state energy when this is not directly observable. However, what can be 
directly measured is the interaction energy Aii^exp between the two electrons in the ground - 
state. This was done for six highly ionized elements, ranging from germanium ( Zex = 32 ) up 
to bismuth (zex = 83), see ref.[l2]. On the other hand, this interaction energy of the two 
ground -state electrons arises within the framework of RST as the difference Ai^^^grp of the 
ground -state energy i?T ^ 13.441 ) and the double value (2M^,c^ ) of the single - particle energy 
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eigenvalues f 13. 251 ). i.e. 

AE^^l.^ = Et- '2M,c^ (3.45) 
= 2 (^Mo - M,) c^ + hcj dVAo(r) ■ ko{r) . 

The comparison of this RST prediction AiJ^grp with the experimental values AEexp displays 
some very instructive features of both the electric and the magnetic interactions ( see table I ). 



Element ( Zex ) 


ASexp [eV] 


a41t [eV] 




f ^ 


e* 


Ge ( 32 ) 


562,5±1,6 


553,0 


1,7 


0,297 


16,8 


Xe ( 54 ) 


1027,2±3,5 


974,3 


5,1 


0,295 


16,58 


Dy (66) 


1341,6±4,3 


1232 


8,2 


0,294 


16,36 


W (74) 


1568±15 


1423 


9,3 


0,247 


16,18 


Bi ( 83 ) 


1876±14 


1661 


11,5 


0,223 


15,92 



TABLE I: Comparison of experimental values AE'exp ( second colmnn ) jlj] with the RST predic- 
tions AE'^grp f 13. 451 ) (third column) for the ground -state interaction energy of helium -like ions. 
The last two columns display the geometric factor f l3.58l ) for the magnetic interactions and the 
reference energy ( 13.511 ') for the electric interactions. Both limit values /q = 0,4 f 16. 431 1 and 
e* 17eV f l3.52l ) for small values of ZexOs(*^ 1) appear to be consistent with the numerical data. 



For an intuitive interpretation of the results of table I it is important to first observe 
that the relative derivation A of experimental data AE^xp and electrostatic RST predictions 
AE'pgrp (fourth column) increases from 1,7% up to 11,5% when the nuclear charge z^-^ ranges 
from Zex = 32 (germanium) up to Zex = 83 (bismuth). This effect can easily be understood 
in terms of the different relativistic behaviour of the electric and magnetic interaction en- 
ergies. Considering first the electric type, one should recall that the single - particle energy 
i?T ( 12.671 ) equals the mass eigenvalue M*c^ ( 13.251 ) which splits up into the matter energy 
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Ed f l2.46al ) and external interaction energy Egg f l2.46cl ) according to 



M,c^ = Mcyi- (zex«s)' = + ^es (3.46a) 

Ed = I (3.46b) 

1 - {zexasf 

Ses = - Mc' = , ^ ^es , (3.46c) 

\2 



see ref. |2]J. Here the non - relativistic limit E^s of the external interaction energy in the 
Coulomb field (^")v4o(r) (EH) is given by 

= < V^o|^'^"^^o|V^o> (3.47) 
= ricj dV^*(^1^^(^1 = -{z,xasfMc' , 
where the non -relativistic ground -state function ipo{r) of the single - particle Schrodinger 

o 

problem f 13. 301 ) coincides of course with the radial wave function -R+(r) f 13. 281 ) 



1 ° , . 11 /z, 



V47r V ^ V«B 



(3.48) 

Ob 

From this result for the external interaction one may conclude that also the internal elec- 
trostatic interaction energy E^j^^ f l2.58al ) of both ground - state electrons is of the following 
form 

Eg) = ^ ^ Ei^^ - AE^^ (3.49) 

y 1 - {ZexOisf 

where E^ is again the non -relativistic limit (^ex^s — > ) of the internal electrostatic inter- 
action energy E^^\ This non -relativistic limit may be simply determined by considering the 
electrostatic interaction energy of the two charge clouds due to the single - particle ground - 
state ipo{r) ( 13.481 ) of the ordinary Schrodinger problem ( 13.301 ): 

^« = //dVdV'M|_^ (3.50) 
= « 17,00725 -2., [eV] . 

o CIb 



Combining this result with the former hypothesis f 13. 491 ) for the electrostatic interaction 

^RST 



energy AE'ifgrp says that the "reference energy" e^, defined through 



1 - {zexasf 

— ■ a41t (3.51) 



Zf,. 
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should adopt the value of ( roughly )17eV for all values of nuclear charge Zex- Now it is just 
this estimate which is confirmed by our numerical calculations ( table I ) within an error of 
a few percent when the elements from z^^ = 32 ( germanium ) up to z^x = 83 ( bismuth ) are 
considered. Indeed, the value Eq to be obtained by a combination of both equations f 13. 491 ) 
and (Einni) 

5 

£0 = = 17,00725... [eV] (3.52) 

8 CtB 

appears as the upper limit for our RST calculations of f 13. 511 ). see table I if the nuclear 
charge z^x is adopted to become smaller ( up to Zex = 2 for the neutral helium; see a separate 
paper ). 

A similar estimate enlightens also the mechanism of the magnetic interactions, whose 
energy - momentum content is characterized by E^^^ f l2.59an . In order to get the magnitude 
of the internal magnetic fields Ha{r), we observe that when a point charge (emitting an 
electrostatic field E at rest ) is moving with velocity v, then there arises a magnetic field of 
magnitude 



H 



E 



(3.53) 



For an extended charge, one would include here some geometric factor being characteristic 
for the special charge distribution, and thus one would arrive at the following relationship 
between the magnetic and electric field strengths 



H 



(3.54) 



Now the kinetic energy part E^ f l3.46bn of the single- particle eigenvalue M^c^ ( I3.46al ) 
suggests to take for the particle velocity in the ground- state 



V 

- ~ Zexas , (3.55) 
c 



and this yields the following relation between electric and magnetic energies AE'^grp ^ Ej^ 
(IHHi) and A^"*"^ 

a4-^) - /.^-%^ • a41t • (3.56) 

Here it is very natural to assume that the magnetic interaction energy IS responsi- 

ble for the discrepancy of the experimental value Ai^exp and the electrostatic RST prediction 
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Ai^^grp, i.e. we rewrite equation f 13. 561 ) as 



Ae!^^^ = AEoxp - AE^'I^T = /* 7-^-7^-^^— 2 ■ ^-^RST ■ (3-57) 

1 - (^exas) 



This then leads us to the expectation that the geometric factor 



2 



1 - (^cx«s)' _ AEexp - AE^'^ 



will depend only very weakly , if at all, upon the coupling constant (zexas)- And indeed, if 
we insert the experimental values Ai^exp and our electrostatic RST results Ai^^grp into the 
right-hand side of equation f 13. 581 ). one just finds the expected weak dependence of see 
for this table I. 

Finally, combining both estimates for the electrostatic part Ai^^grp ( 13.511 ) and the mag- 



netic part Ae!^^^ (ESni), we arrive at the general form of the interaction energy AE, 



cxp 



AEexp = AE^^l^ + A^S (3.59) 

\2 



1 - (ZcxttS 



|2 1 - (2;exas) 



This general result contains two slowly varying functions of the coupling constant (-ZexCts), 
which adopt their maximum values for small nuclear charge numbers [z^^ 2 ) . The 
maximum value of £0 of the reference energy has already been determined in equation 
f 13. 521 ). but for determining the maximal value (/q, say) for the geometric factor 

/o=0,4 (3.60) 

one has to go deeply into the details of the magnetic interaction mechanism ( see the dedu- 
cation of equation f 16. 431 ) below). As we shall readily see, the crucial point here is the fact 
that the magnetic exchange forces do not vanish since there exists a non- trivial exchange 
current h{r), producing the magnetic exchange effects such as the magnetic exchange energy 
E^^ f l2.59bn . whereas for the ground - state there is no electric exchange energy E^^ f l2.58bl ) 
on behalf of the vanishing of the scalar exchange potential BQ{f)\ 

IV. MAGNETIC CORRECTIONS 

The preceding estimate of the magnetic energy contributions hinted strongly upon a 
magnetic explanation for the discrepancy between the experimental values Ai^gxp and the 
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electrostatic RST predictions Ai?^grp (table I). Therefore it appears now mandatory to 
explicitly compute the magnetic corrections, thought to be responsible for the observed dis- 
crepancy. Through this procedure one will then obtain also more accurate RST predictions. 
In the course of such an elaboration of the magnetic corrections it should also become clear 
in which way the upper limit /g ( I3.(i()l ) of the slowly varying function ( 13.581 ) comes 
about. Furthermore, the role of the magnetic exchange effects , being induced by the "mag- 
netic" exchange potential B{r), has to be clarified. Indeed, it will readily turn out that this 
"magnetic" exchange potential B{r) cannot vanish for the reason of ground - state symmetry 
( i.e. isotropy ), which is in sharp contrast to the missing of the "electric" exchange potential 
Bo{f). Actually, the "magnetic" exchange energy due to B{r) will turn out to be twice the 
magnetostatic contribution due to Aa{r), which yields a higher precision of the RST pre- 
dictions by roughly one order of magnitude (compare the RST results of tables I and II). 
Consequently, the experimentally supported inclusion of the "magnetic" exchange effects 
due to B{f) must be viewn as a strong confirmation of the non- ahelian RST construction. 
Thus it becomes now necessary to work out the field theory of atomic magnetism in great 
detail. 

The desired magnetic corrections do emerge in two places: 

(i) as a small change M^^^ of the mass eig envalues M„ (E3i) - (ESEl) , 

(ii) as an additional contribution E^^^^ f 12. 621 ) of the interaction energy Eq f 12. 601 ). 
Therefore the magnetic perturbation scheme must consist in first expressing the mass 
changes Mi'"^^ and energy change E^^^"* in terms of the wave functions ipai'r) as solutions 
of the relativistic eigenvalue problem ( I3.4al ) - ( I3.4bl ). and then substituting for these solu- 
tions their electrostatic approximations ^/'a(^) as solutions of the truncated system ( I3.6al ) - 
f l3.6bn . In other words, one computes the value of the total energy functional i^x upon the 
approximative solutions ^/'a(r) place of the exact solutions V'a(r) original system 

(mi) -(Eli). 

A. Magnetic Mass Corrections 

The point of departure for obtaining the "magnetic" mass corrections is the original exact 
form ( I3.4al ) - ( I3.4bl ) of the two - particle problem. The corresponding eigenvalues Ma are 
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splitted into their electrostatic approximations Ma and their magnetic corrections mIT'^^ 

Ma = Ma + Mjr^\ (4.1) 

Quite generally, the exact mass eigenvalues Ma can be expressed in terms of the correspond- 
ing eigenfunctions ipair) by multiplying both sides of the eigenvalue equations ( I3.4al ) - ( I3.4bl ) 
by V'a(r) ^^"^ integrating, observing also the normalization relations f 13. 351 ). The mag- 
netic mass corrections will then be due to those terms containing the vector potentials 
Aair) = {'■"•'^•'(r)} and B{f) = {B^(r)}, which have been omitted for the electrostatic ap- 
proximation. However, for discussing these magnetic terms it is instructive to reunite them 
with the electrostatic interaction terms and thus to consider the interelectronic interactions 
as a whole. This admits us to resort to the covariant Maxwell equations for carrying out 
the magnetic perturbation procedure where afterwards the separation of the electric and 
magnetic effects for obtaining the desired result Ma"^^^ will present no problem. 

In this sense, the united mass corrections AM^ may be deduced from the eigenvalue 
problem ( I3.4al ) - ( I3.4bl ) in the following form: 

zi ■ AMic^ = -he J d¥{(i)A;^(f) ■ ^^W{r) + h^{r) ■ E'^(r)} (4.2a) 

Z2 ■ AMac^ = -he j d¥{(%(f) ■ (^U'^(f) + h*^{r)-B^'*{r)] . (4.2b) 

Obviously the mass corrections are built up by an electromagnetic part m'^'^'' and an ex- 
change contribution Ma'^^^ 

AMa = M^""") + Mf ^) , (4.3) 
with the self-evident arrangements 

zi ■ Mi'^^^c" = -he [ dV ^%{r} ■ ('U^(f) (4.4a) 



Z2 ■ M^^'^^c^ = -he J d¥ (%(f) ■ (lU'^(f) (4.4b) 
zi ■ M^^^^e^ = -he J dV h^{f) ■ B^ir) (4.4c) 
Z2 ■ Mf ^^c^ = -he J dV h*^{f) ■ B^'*{r) . (4.4d) 

Here, the electromagnetic mass corrections Mjf'^'^ can be further split up according to 
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with the electric part Ma*^^ being given by 

. m['^c^ = -he j d¥ ^%o{f) ■ (2Uo(f) (4.6a) 

^2 ■ M2^'^c2 = -he j d¥ ^%o{f) ■ (^)Ao(f) , (4.6b) 

and similarly for the magnetic part iVfi'"^ 

zi ■ M^^^c^ = -he j d¥ ^%ir) ■ (^U^ (f) = j d¥ A^i(f) ■ A2{f) (4.7a) 

Z2 ■ M!p^e^ = -he j d¥ ^%j{r) ■ (^U-' (f) = J d¥ hir) ■ Mr) . (4.7b) 

In a similar way, the exchange corrections Mji^^^ f l4.4cl ) - (0331) can also be subdivided 
into their "electric" parts Ma^^ and "magnetic" parts Ma^^ 

Mihg) = jvff ) + M^») , (4.8) 

i.e. we put 

zi ■ Mf ^c^ = -he J d¥ ho{r) ■ B^f) (4.9a) 

Z2 ■ Mf ^c^ = -he j d¥ h*ir) ■ B^*{r) (4.9b) 

zi ■ M^{'^e^ = -he j d¥ hj{r) ■ B^{r) = j d¥ h{r) ■ B{r) (4.9c) 

. M^^^e^ = -he j d¥ h*{r) ■ B^*{r) = ^ j h*{r) ■ B*{r) . (4.9d) 

The physical meaning of these mass corrections becomes now evident when they are 
written in terms of the curvature components -F^^^, and Gf^^. Such a reformulation of the 
mass corrections can easily be attained by eliminating the currents and /i^ from the 
original definitions ( I4.4al ) - ( I4.4dl ) in favour of the curvature components, namely just by 
means of the covariant Maxwell equations ( I2.36al ) - ( I2.36dl ). Now, at this point we resort 
to an additional approximation which considerably simplifies our magnetic perturbation 
approach, namely the linearization of just that gauge field dynamics ( I2.36al ) - ( I2.36dl ). 
This linearization means concretely that 

(i) the cuvature components -F"^^^, Gf^u ( I2.31al ) - ( 12. 3 lei ) become truncated into 

= d,A\ - d,A\ (4.10a) 
G^, = d^B,-d,B^ (4.10b) 



32 



by simply omitting the non- linear terms due to the non-abelian character of the 
original gauge group U{2); 

and 

(ii) the non-abelian Maxwell equations f l2.36an - f l2.36dn become similarly truncated to 
their linear form 

d^T"^^ = -Anas k\ (4.11a) 

d>'G^, = A7TasK. (4.11b) 

As a consequence of these simplifying assumptions, one finds for those four- vector prod- 
ucts dertermining the mass corrections Mi""^ and iVff (WJ^) - (\4Ad\): 

^%{f) ■ (2U^(r) = -J—d^\WF,,{r) ■ (2U-(f)| + J—mF,,{f) ■ (2)F'^^(f) (4.12a) 
(%(f) ■ (lU'^(f) = -^9''|(2)F^,(f) ■ (iU'^(f)| + -^(2)^^,(f) ■ (i)F'^'^(f) (4.12b) 
h,ir) . B^{r) = -^d>^{G%{r-) ■ 5'^(f)} - J—G*^^-^ (f) ■ G,.(r) . (4.12c) 

When this is substituted back into the correction formulae ( I4.4al ) - ( I4.4dl ). they appear in 
a new form by means of Gauss' integral theorem, namely as 

. "^c^ = Z2 ■ M^'"'^^ = — [ d¥ (^)F„,(f) ■ (^)F^-{f) (4.13a) 

8nas J 

. ^^c^ = Z2 ■ M^^'^c^ = / d¥ G*''" (f) • G.,{r) . (4.13b) 

Svras J 

This result explicitly demonstrates that the mass corrections AM„ ( KM ) - ( K2Bi ) for both 
particles ( a = 1, 2 )are actually identical ( zi ■ AMi = zi ■ AM2 ). 

But once that "Lorentz invariant" form ( I4.13al ) - ( I4.13bl ) of the mass corrections is 
known, it is self- suggesting to split them up into two contributions, which are offered by 
themselves through the introduction of the electrostatic and magnetostatic field strengths 
ESf) = {^"^E^ir)} and Ha{r) = {^"Wif)}, i.e. we put 

WE^(f) = ('')Foj(r) = -dj'^^%{f) (4.14a) 
Wifi(f) = ^e'\ ^"^F^'ir) = 9j'^U'(f) , (4.14b) 
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or in three -vector notation 

Ea{^ = -V("Uo(f) (4.15a) 
Hai^ = V X ^^(r) . (4.15b) 

Now by this arrangement, the Lorentz invariant product of curvature components deter- 
mining the electromagnetic mass corrections Mjf™'^ f l4.13an reads in three -vector notation 

(i)F^,(r) ■ (2)F'^-(r-0 = 2 [#i(f) ■ #2(r) - E^{f) ■ E2(f)] , (4.16) 

and clearly this yields a natural splitting of those electromagnetic mass corrections f l4.13al ). 
namely just the former equation f l4.5l ) with the following identifications for the electric part 

■ m['^c^ = Z2 ■ M^'^c^ = -p— [ d¥ Ei{r) ■E2{f) , (4.17) 

47ras J 

and similarly for the magnetic part 

5i • M^f'^c^ = Z2 ■ M^r'^c^ = — [ d¥ #i(f) ■ #2(rO . (4.18) 

4vras J 

Observe here the curious fact that the magnetic part ( 14.181 ) enters the electromagnetic mass 
correction Mjf"^^ ( \A.b\ ) with the opposite sign in comparison to the electric part ( 14.171 ) which 
is a consequence of the Lorentz invariance of the product of the field strenghts f 14. 160 ! 

Clearly, it is self- suggesting now to treat the exchange corrections m!i^^^ f l4.13bn in 
a quite similar way. This means that one introduces an "electric" exchange vector field 
X{r) = {X^{r)} (as the exchange counterpart of E{r) ( I4.14al ) )and also a "magnetic" 
exchange field Y{f) = {Y^{f)} (as the exchange counterpart of H{f) ( I4.14bl ) )through 

X^{f) = Gojif) = -d,B,{f) - —B,{r) (4.19a) 

Am 

= l^^'^Gkii^ = e'\dkB\r) , (4.19b) 
where the exchange length parameter om is given by 

In three -vector notation, the relations ( I4.19al ) - ( I4.19bl ) read 

X(f) = -V5o(f) + —B{f) (4.21a) 
Cm 

Y{r) = V X B{r) , (4.21b) 
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and thus obviously represent the exchange analogue of the corresponding electromagnetic 
relations f l4.15al ) - f l4.15bn , however with the difference that there emerges now a typical 
length parameter om f l4.2UI ) which gives an inherent measure for the spatial range of the 
exchange effects. In terms of these new exchange fields, the Lorentz scalar for the exchange 
mass corrections Mj^^^^ f l4.13bl ) reads 

G*^"" (f) ■ G^^(f) = 2 \Y{f)* ■ Y{r) - X{r)* ■ X(f)j , (4.22) 

which is of course again the exchange analogue of the corresponding electromagnetic relation 
( 14.161 ). As a result, the exchange mass corrections M^^^^ (EIM) can ultimately be written 

as 

zi ■ Mf ^^c' = Z2 ■ Mf ^^c' = f d¥ [y (f)* ■ F(f) - X{r)* ■ X{r}] . (4.23) 

47ras J L J 

Though this result looks formally quite analogous to the electromagnetic case, one never- 
theless cannot relate here the "electric" exchange masses Mi''^ (EM) - (USE) to the "elec- 
tric" exchange vector X{f) and analogously the "magnetic" corrections Mjf^ f l4.9cl ) - f l4.9dl ) 
to the "magnetic" exchange vector Y{f) as it was done for the electromagnetic situation 
( 14.171 ) - ( 14.181 ). Instead, the correct relationships for the "electric" exchange mass must 
look as follows: 

zi-Mf^c^ = i2-Mf ^c^ = — [ d¥ X(r)*-X(r) + -^^ [ d¥ B(f)*-B(r) . (4.24) 

47ras J 47ras au J 

Thus both "electric" exchange corrections are again identical, as it is the case also for the 
electromagnetic subsystem f l4.17l ). However, they do not coincide with the "electric" part 
due to X{r) of the total exchange corrections Mi^^"* ( 14.231 ). 

A similar effect occurs also with the "magnetic" parts Ad'if'' of the exchange mass cor- 
rections M^^"^ f 14. 231 ). Starting from the original definitions ( I4.9cl ) - ( I4.9dl ) one finds again 
that both "magnetic" contributions are identical 

zi ■ Mi'^c^ = Z2 ■ M^'^c^ = f dV Y{r}* ■ Y{r) - f d¥5*(f) ■ B{r} , (4.25) 

4vras J 47ras J 

but they do again not coincide with the "magnetic" part ( due to Y{f) )of the total exchange 

masses Mi^^^ (EM)- However, if both the "electric" and the ' 'magnetic" , contributions 

( 14.241 ) and ( 14.251 ) are added up separately for either particle ( a = 1, 2 ), one just finds the 

total exchange corrections Ma^^^ ( 14.231 ). Thus the electromagnetic and exchange subsystems 

are found to differ in the way in which the total corrections are distributed upon the "electric" 

and the "magnetic" subsystems. 
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B. Linearized Gauge Field Equations 



Obviously, the introduction of the three - vector notation enables one to separate uniquely 
the electric corrections ( being already included in the electrostatic approximation ) from the 
magnetic corrections which will be treated subsequently by use of an adequate perturba- 
tion approach. Consequently, one would like to write down now the gauge field equations 
f l2.36an - f l2.36dn in three - vector notation where one simultaneously restricts oneself to their 
linearized form f l4.11al ) - f l4.11bn . This implies that the magnetic corrections will be taken 
into account only in their lowest -order approximation. However, the consistent lineariza- 
tion of the gauge field equations represents a certain problem which requires now an extra 
discussion. 

First consider the electromagnetic subsystem, for which one concludes from the linearized 
Maxwell equations f l4.11an that the three - current ka{r) = | '^^^A;-' (r) } must have vanishing 
divergence ( a = 1, 2 ) 

V ■ kair) = . (4.26) 

Thus the linearization of the (non-Abelian) Maxwell equations f l2.36al ) - f l2.36bn must im- 
ply the neglection of the right-hand sides of the source equations f lTTO l - f lCTbl V Fur- 
thermore, when the curvature components -F^^j, are expressed through the connection com- 
ponents A"-^ f l4.10al ). that equation f l4.11al ) yields the corresponding Poisson equations 

A(''Uo(f) = 47ras ■ (%(f) (4.27a) 
AL(f) =47r«s-4(r) , (4.27b) 

provided one subjects the magnetic vector potentials Aa{r) ( a = 1,2 )to the Coulomb gauge 
condition 

V • X(f) = . (4.28) 

Observing here the usual boundary conditions at infinity (r — > oo), the solutions of the 
Poisson equations f l4.27al ) - f l4.27bn are easily found as 

(«)Ao(f) = -as / d\-" ^.^ML^ (4.29a) 
J \f— r'\ 

A,{f) = -as [ d¥' M!) , (4.29b) 
J \r — r'\ 
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where the spatial part f l4.29bl ) is easily seen to meet with the gauge condition f 14. 280 . 
Thus it is evident that the linearized electromagnetic subsystem is governed by the usual 
( Abelian )Maxwellian structure, which may be expressed also through the field equations 
for the field strenghts Ea{r), Ha{r) (Eni) - (ElEl) : 

V ■ E^(f) = -47ras ■ ^%(^) (4.30a) 

V ■ Ha{r) = (4.30b) 

V X Ea{r) = (4.30c) 

V X i?,(r) = -47ras • 4(r) . (4.30d) 

However, the exchange subsystem has a somewhat different structure. This became 
already obvious through the introduction of the "electric" exchange field X{r) f l4.21al ) 
which is not simply the gradient of the "electric" exchange potential Bo{r) but contains 
also the "magnetic" exchange potential B{f)\ Nevertheless, one finds from the linearized 
equations ( I4.11bl ) again the ordinary Poisson equation for -Bo(^) 

A5o(f) = -47r«g/i*(f) (4.31a) 
Boirl =asj dV'^^^ , (4.31b) 

which is thus revealed as not being affected by taking into account the magnetic corrections 
( similar as for the electric potentials ^""^AQ^r) f l4.27an ). Observe also, that for the deduction 
of the exchange Poisson equation f l4.31al ) we imposed a Coulomb - like gauge condition upon 
the "magnetic" exchange potential B{f) 

V ■ B{f) = , (4.32) 

in close analogy to the magnetic counterpart Aa{r) f 14. 281 ). Strictly speaking, such gauge 
fixing conditions as ( 14.281 ) must not be imposed upon the fields B^ because they lost their ge- 
ometric status as gauge potentials in the pro cess of Abelian symmetry breaking and therefore 
obey a homogeneous transformation law j2l|. However, it is easy to see that the requirement 
f 14. 321 ) is consistent with all the other static field equations for the exchange variables Bo{r), 
B{r). 

It is true, the "electric" part of the exchange field equations f 14. 31 al l - f l4.31bn is not too 
much different from its electric counterpart ( I4.27al )and ( I4.29al ): however, for the exchange 
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vector potential B(r) one finds from tlie linear equations f l4.11bn 

AB{f) + 4- ■ B{r) = -47rash*{r) - —VBo{r) . (4.33) 

As a consistency test of this equation one forms here the divergence of the left- and right- 
hand side and arrives just at the Poisson equation ( I4.31al ) for the exchange potential Bq^v), 
provided the divergence relation for B{r) f 14.321 ) is respected together with the following 
source equation for the exchange current h{r): 

^ .h*{f^ = J_hl{r) . (4.34) 



However, this is nothing else than the continuity equation for the exchange four -current h 



d^'hlir) = (4.35) 

which itself is a consequence of the linear exchange field equations f l4.11bn and, on the 
other hand, coincides with the former source equations f l2.4(Jcl ) - f l2.40dn if the right - hand 
sides vanish (to be justified below). In this way, one actually attains a consistent linear 
approximation of the non-Abelian (and therefore non- linear) gauge field equations. 

Of course, the solutions B{r) of the exchange field equation ( 14.331 ) must have a somewhat 
other shape than their magnetic counterparts Aa{r) ( I4.29bl ): 

B{f} = as [ d¥'— ^ ■ cos (t^ll] . Sh*{f') + V'5o(f ')) . (4.36) 

J \r-r'\ \ au J y 47rasaM J 

Thus the effect of the exchange length is to suppress the magnetic exchange vector 
potential B{r) if the exchange current h{r) is smoothly spread over a spatial domain much 
larger than the exchange length om- But if h{r) is well- localized within such an "exchange 
domain", the vector potential B{r) is non -zero, but fades away as ■ cos On 
the other hand, for ayi ^ oo the "magnetic" exchange equation f 14. 331 ) degenerates to an 
ordinary Poisson equation of the kind ( 14.3 lal ) or ( I4.27al ) - ( I4.27bl ) with the corresponding 
behaviour of the solutions. Thus the "magnetic" exchange effects are found to be of a rather 
different type which may again be expressed more concisely by the following "exchange 
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equations" : 

V ■X{r) = A7rash*{r} (4.37a) 

V X X{r) = —Y{r) (4.37b) 

V ■ Y{f) = (4.37c) 

V X F(f) = 47ras ■h*{r) + ^- B{r) + — • V5o(r) , (4.37d) 

which adopt the ordinary Maxwelhan form f l4.30al ) - f l4.30dn for infinite exchange length 
(om — > oo). As a consistency test for the hnearization procedure, apply the divergence 
operation to the last equation ( I4.37dl ) and find just the source equation for the exchange 
current h{r) ( 14.341 ) by use of the Poisson equation ( I4.31al ). 

C. Mass Corrections and Densities 

Once both the electromagnetic potentials |'''*Mo(r*) ; 74a(r)| and the exchange potentials 
|i?o(^) ; -B(r)| are known in terms of the charge and current densities, it becomes possible 
to eliminate these potentials completely from the mass corrections and to express the latter 
objects exclusively in terms of these physical densities. Clearly, this then represents a 
considerable technical simplification because one can express the desired mass corrections 
directly in terms of the wave functions whose link to the current densities is well-known, 
see equations ( 13. Sal ) - ( I3.8bl ). Thus to begin with, reconsider the electric corrections M'f^ 
f 14. bah - f l4.6bn and substitute therein for the electrostatic potentials ^"''^AQ{f) their general 
form f l4.29al ) in order to arrive at 

Of course, the charge densities ("^fco(r) here be further expressed in terms of the Dirac 
wave functions ipai'f^) according to the former relations f l3.8al ). Observe also, that the 
present result of the electrostatic interaction energy Za ■ Mjf^c^ now admits two different 
interpretations, namely either as an instantaneous Coulomb interaction between the two 
charge clouds ^"•^ko{r) or as the interaction energy {not self- energy! ) of the electric field 
modes Ea{r) emitted by the charge clouds, see equation ( 14.171 ). This bilinear (instead of 
quadratic ) construction for the electromagnetic interaction energy of the particles has a cer- 
tain tradition in the literature and emerges also quite naturally in the RST formalism 
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A similar effect is obtained also for the magnetic corrections Mi*"-* f l4.7al ) - f l4.7bn which 
by means of the vector potentials Aa{r) f l4.29bn may be recast into the following form 

z, ■ Mi^^c' = z, ■ Mt^c' = -e' [[ d'rd'f'^^^^r^^p- . (4.39) 



Thus the magnetostatic interaction energy Za-Mi^'^c^ can also be interpreted as being due to 
either an instantaneous direct interaction of the currents ka{r) or as the interaction energy 
of the magnetic field modes Ha{r) f 14. 181 ). It is also interesting to remark, that the the sum 
of the electric and magnetic corrections appears as the integral of a four -vector product, 
namely by adding both equations f 14. 381 ) and f l4.3 



za ■ Mi'^c" + za ■ Mjr^c^ = z^ ■ Mjf"'^ (4.40) 

= ff d'fd^'f' 



.2 II ,3^,3j''\i^-^''^k'^ir') 



{a = 1 or 2) . 

Now it seems natural to expect that the exchange field system will display an analogous 
structure. This however is true only for its "electric" component, but not for the "magnetic" 
one. In order to see this more clearly, first recall that the "electric" exchange potential BqI-t) 
(KSM) - (EMH) is of the same structure as its electric counterpart ^"'^Ao(r) ( I4.29al ) and 
therefore yields for the "electric" exchange mass Ma^^ ( I4.9al ) - ( I4.9bl ) 

z, . Mi'^c' = z, ■ M!t^c' = -e' [[ rfV-rfV-M^^^P . (4.41) 



Indeed, this looks very similar to the electrostatic corrections Mi"^ (1123) where the charge 
densities ^'^^ko{r) play the part of the exchange densities ho(r), K^{r) and a change of sign 
does occur additionally. This lowering of the electrostatic interaction energy ( 14.381 ) by the 
exchange energy f 14. 411 ) is due to the fact that the present RST formalism is equivalent to 
the antisymmetrized product states of the conventional Hartree - Fock approach. 

However, the "magnetic" exchange system appears in a somewhat different shape which 
traces back to the modified Poisson equation ( 14.331 ) for the "magnetic" exchange potential 
B[r). More concretly, introducing the solution ( I4.3(il ) for B{r) into the exchange correction 
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formulae f lQcl ) - f l4.9dn yields 



cos(^U(f) ■K*(r") 
h ■ m['^c^ =e^ £f£f' ^ ^, (4.42a) 



,hc rr cos f ^U(f) . v'So(r 



47ra 



M 



COS f^^) /1(f) ■K*(r") 
Z2 ■ M^^^c^ = / / dV-d^f' ^ ""7 ^, (4.42b) 



47raM 

Obviously, the presence of the terms containing the gradient of the "electric" exchange po- 
tential Bo{r) does not allow us to recast this result again into the form of an integral over 
the four -vector product h*^{r) ■ h'^{r), as it was the case for the electromagnetic analogue 

(loni). 

After the energy contributions of the electric and magnetic type have been discussed 
in detail, one can now render more precise the meaning of the "magnetic corrections". 
Obviously, the purely electric corrections Mjf'^^ (to be included into the electrostatic ap- 
proximation ) are given by all the mass corrections of the "electric" type ( 14.381 ) and ( 14.411 ) 

whereas the magnetic corrections ( to be superimposed as perturbations over the results of 
the electrostatic approximation ) are given by the correction terms of the "magnetic" type 
(KM) and (|432i) - ( \4A2R ) : 

The lowest -order approximation of such a perturbation scheme consists then in inserting the 
solutions ipai^) of the electrostatic approximation system f l3.6al ) - f lS.Gbl ) into these magnetic 
correction functionals ( 14.441 ). 
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D. Magnetic Energy of the Gauge Field 



Besides the magnetic contributions M^^^ as part of the mass eigenvalues M^, there 
occur further magnetic constituents of the total energy ii^T? namely those being due to the 
gauge field energy Eq f 12. 531 ). Since the gauge field energy Eq f 12. 531 ) itself splits up into 
the electromagnetic part -Er f 12. 541 ) and exchange part Eq f 12. 551 ). one may separate again 
both parts into the electric and magnetic type f l2.57al ) - f l2.57bn . Here the electric part E^^ 
f l2.58an is naturally included in the electrostatic approximation and is found now to agree 
just with the electrostatic mass corrections m'^'' f l4.17l ): 

= Za ■ Mi'^c^ , ( a = 1 or 2 ). (4.45) 

This is the reason why the subtraction of Eq^^ from the sum of mass eigenvalues rescinds the 
double counting of the electrostatic energy E^^'^ for the total energy Ej' f 13. 391 ). Similarly, 
the magnetic part E^^^ of the electromagnetic gauge field energy Eji f l2.57al ) reads in terms 
of the magnetic field strenghts Ha{r) 

Ef) = J dViJi(f) ■ #2(r) (4.46) 

and thus is the negative of the magnetostatic mass corrections Mi"^ ( KW ) 

= -z^ . Ml^^c^ ( a = 1 or 2 ) , (4.47) 

in contrast to the corresponding electric counterpart ( 14.451 ) . Therefore the addition of Eq^^^ 
to the sum of mass eigenvalues rescinds the double counting of the magnetostatic energy 
E^^^ for the total energy (E3i)- 

In an analoguous way, the exchange energy Eq f l2.57bl ) may also be split up into the 
contributions of "electric" and "mag netic" type according to eJ!"^ r l238Hl ) and E^^'^ r l239El ). 
The "electric" part i?^^ however cannot be identified here with the mass corrections m2^^ 
f l4.9al ) - f l4.9bl ). as it was possible for the electrostatic approximation, but the exchange 
vector potential B{r) has to be retained so that one first finds 

eP = — ^ I d^fB*(f) ■ Bir) + f dVV5*(f) • VB*Jr) . (4.48) 

47rasaM J ^vras J 

Clearly, when the exchange vector potential B{r) is neglected here, one returns to the 
truncated form due to the electrostatic approximation. On the other hand, eliminating the 
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"electric" exchange potential Bo{f) from the present result for f 14.481 ) in favour of X(r) 
by means of its definition f l4.21an yields now the generalization of the former electrostatic 
approximation, namely for both a = 1 or a = 2: 

^(h) ^ _hc_ I d3^^*(^) (449) 
47ras J 



Za ■ Mi'^^c' + — ^ / d¥5*(f) • Bir) . 



Thus it is seen that the "electric" exchange energy E^^ receives an additional contribution 
in comparison to its mass correction counterpart M^'^'' f l4.9an - f l4.9bn when the "magnetic" 
exchange potential B{f) is not zero. This additional contribution (i.e. the integral on 
the right - hand side of f 14. 491 ) ) must be taken into account for building up the "magnetic" 
energy correction AEq^^^ below; whereas the first contribution, i.e. the "electric" exchange 
mass Ma^"^ ( 14.411 ). is already contained in the electrostatic approximation of Et ( 13.421 ) and 



must therefore be omitted for the magnetic part of AEi™^^ 



In a quite analogous manner one can treat the "magnetic" exchange energy 
which is nothing else than the exchange analogue of the magnetostatic field energy -E^^ 
( 14.461 ). Comparing this again to the "magnetic" exchange masses Mi^^ (11211) yields the 
relationship 

Ei^^ = Za ■ M[3)^ + I d¥5*(f) ■ B{r) . (4.50) 

4:Tcasaf^ J 

Here it is interesting to observe that the energy corrections of the electric type, namely E^^^ 
f 14. 451 ) and Eq^ f 14. 491 ). contain the corresponding mass corrections Mjf^ and Ma'^^ with 
opposite sign so that they actually become subtracted from the sum of mass eigenvalues Ma 
in order to build up the total energy E^ f 13. 391 ). As a result, the double counting of these 
corrections of the "electric" type is actually avoided, as explained in connection with the 
electrostatic approximation (see the discussion below equation ( 13.411 ). 

However, in contrast to this, the corrections of the mag netic type, i.e. E^^^ ( WJ7\ ) and 
E^f^ ( WW\ ) both contain the " magnetic" mass corrections Mi™'* and Mjf^ with the "false" 
sign; but since E^(^^^ (121121) has to be added to the same terms occuring in the sum of 
mass corrections Ma for building up the total energy £'t f 13. 391 ) the double counting of the 
"magnetic" energy contributions is rescinded in just the same way as it is the case with the 
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"electric" contributions. For the magnetic energy correction AE^^^ we therefore find 

AEg"^) = - - / d¥5*(f) ■ B{r) (4.51) 

Ana 



= -4 ■ Mi-)c2 - 4 ■ Mi^)c^ - — ^ / d'fB*{r) ■ B{r) . 

Observe here that half of the integral on the right - hand side emerges for both the "electric" 
part E^^'^ ( 14.491 ) and the magnetic part E'^^ ( 14.501 ) and thus must enter the energy correction 
AEjf (HSU) with its double value. 

Summarizing, the total energy correction AE'^^^ of the "magnetic" type is now found 
as the sum of all those terms in the total energy E'^ ( 13.391 ) which were omitted for the 
electrostatic approximation Et ( 13.421 ). i.e. 



A4-g) = J]i,Mi"^^)c2 + AEjr^^ 

a=l 

= \ t - ^ / 'i'fS'if) . B(r^ . (4.52) 

Here the "magnetic" masses are the proper sum of the magnetostatic mass 

(HSni) and its exchange counterpart Mi^^ ( I4.42al ) - ( f4.42bl ). see equation ( lOl ). Clearly, 
the expectation is now that, when the magnetic energy correction AE^^^^ (11321) is added 
to the electrostatic approximation Et f 13. 421 ) one will arrive at a more precise numerical 
prediction for the atomic energy levels E't. However, before this expectation receives its 
validation, it is very instructive to convince oneself of the physical correctness of the RST 
picture of the atomic magnetism by considering the interaction with an external magnetic 
field ( Zeeman effects). 



V. EXTERNAL MAGNETISM 

It may seem somewhat strange that the magnetic contributions Eq^^^ are entering the 
total energy i?T f 13. 391 ) with the opposite sign in comparison to the electric contributions 
AEg''). This is the more amazing as the magnetostatic energy content E^^^ i \2.f)\)ai \ is 
also given by the positive product of the magnetostatic fields Haif), quite analogously to 
the electrostatic counterpart e'^^ ( I2.58al ). However, the formalism of minimal coupling 
f l3.5al ) - f l3.5bl ) together with Lorentz invariance inevitably leads to that minus sign for 
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the magnetostatic corrections as part of the total energy E^; and therefore one may wish 
to have some additional supporting argument that the magnetic interactions are actually 
implemented in the RST dynamics in the right way. Such an additional argument can be 
put forward by considering the interaction of either of the two electrons with an external 
constant magnetic field Hex- Indeed, the covariant derivatives f l3.5al ) - f lS.Sbl ) most clearly 
display the fact that the corresponding vector potential Aex{r) 

V X Aexir) = Hex (5.1) 

acts upon the a-th particle in just the same way as does the vector potential ^'^^A{r) (b ^ a) 
due to the other particle. Therefore, we can test the correctness of the magnetostatic inter- 
electronic interactions by simply inspecting the interaction with an external source emitting 
the constant field Hex ( Zeeman effect, see e.g. js^] ). 

Now it is well well known that the interaction of a bound electronic system with an exter- 
nal magnetic field is phenomenologically described by including into the ( non - relativistic ) 
Hamiltonian H an interaction term Hint of the following form 

Hint = —f^J ■ Hex ■ (5.2) 

Here the operator of the total magnetic moment /Tj of the system is composed additively of 
the orbital and spin parts 

/^j = /^L + /?s , (5.3) 

which themselves are proportional to the corresponding angular momentum operators L and 
S resp., i.e. 

= (5.4a) 

k = -2^1 (5.4b) 

e^h g magneton I . 
^ 2Mc ' ^ J 

Thus, e.g. when the Russell - Saunders coupling occurs, the total magnetic moment fij ( \b.3\ ) 
can be written as 

lij = -93-^J, (5.5) 
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where the gyromagnetic ratio is given by the Lande g- factor. In any case, the interaction 
Hamiltomian Hint ( 15.21 ) imphes the existence of an interaction energy Emt 

Ei„t = < i^\HinM > (5.6) 

which is experimentally confirmed very well by atomic spectroscopy. This says however 
that, when the magnetic interactions are correctly incorperated into the RST dynamics, the 
non- relativistic formula ( I5.fil ) for the external magnetic energy Eij^t must be also deducible 
from our RST results. 

Indeed, such a deduction can easily be attained in the following way: restricting ourselves 
for a moment to a single particle (a = 1) with normalized four -current k^{f) f 13. 261 ). its 
external magnetic energy say) is deduced either directly from the one -particle 

equation (|S3H) with A^^ = = 0, or from the one -particle version of the energy func- 
tional Et {EM)- 

M,c^ + E^") . (5.7) 

Indeed, if one could omit here the external term E^l^^ f 12. 480 (see below), so that the field 
energy i?T equals the mass energy M^c^ j2l|, both methods would yield the same result: 



Mi^^c^ = hc d^fk{r) ■ Aec{r) . (5.8) 



This is exactly the way in which the considered particle would also magnetically interact 
with the other (not considered) particle, cf. ( I4.7al ) - ( I4.7bl ). so that there is actually no 
difference between external and internal magnetism. However the presumed omission of that 
external term E^s^^ in equation ( 15.71 ) is actually justified for the presence of a homogeneous 
external field H^x- The reason for this is that in this case the magnetic volume integral 
in equ. ( 12.481 ) may be converted to an ( ill - defined ) 2 - surface integral at spatial infinity 
(r — i> oo). This surface integral has to be conceived as the (infinite) energy content due 
to the (infinite) external source emitting the homogeneous field Hex', and therefore that 
external term E^s^^ is to be omitted when one considers the energy of a localized particle. 

Now the vector potential A^xir) due to a constant magnetic field H^x is given (apart 

from a gauge transformation) by 

Aexif^ = -\ (r X Hex) , (5.9) 
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and therefore the external magnetic energy f l5.8l ) is found to be of the following form 

Mi"^^c' = -fl,-HitK (5.10) 

Here the magnetic moment /Ij is given by 

1 /".s. 



/ij = --e*y dV(^fx /c(f)j , (5.11) 

where e^, denotes the elementary charge unit and the physical field Hex^^ is related to our 
geometric notation ifex through 

He. = l^^if ) . (5.12) 
Next, recall that the Gordon decomposition of the four -current ^^j^ of a single Dirac 



particle reads 



2l| 



where the objects cr^,y are the Spin(l, 3) generators, i.e. 

(^i.u = ^ bi^lu] ■ (5.14) 

Observe here that we are allowed to resort to the Dirac current f 15. 131 ) of a free particle, 
because the magnetic energy ( 15.101 ) already contains the external field ifex linearly and we 
are satisfied with the deduction of that first -order approximation ( 15.61 ). Obviously, the 
Gordon decomposition ( 15.131 ) splits up the Dirac current ^^^j^ into a sum of two parts, 
namely the drift part *^°'^/c^ 



and the polarization part ^^^A;^ 



% = (^^/^) , (5.16) 



where both parts obey a separate continuity equation, i.e. for the free particle: 

d>'^% = (5.17a) 
= . (5.17b) 
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For the Gordon decomposition of a composite system see ref. [2]|). Now insert the sphtting 

k, = ^% + ^% (5.18) 

into the definition of the total magnetic moment /2j f 15. Ill ) and find its analogous splitting 
as 

N = ^L + f^s, (5.19) 
with the orbital part fii being given by 

/iL = —e^J d¥ (fx fc,(f)) , (5.20) 

and similarly for the spin part /Is 

/Is = -Je, / d¥ ff X kp{f)) (5.21) 



2 

fcd(r) = {('^V},4(f) = {(*')A;^} 

But clearly, the splitting of the total magnetic moment /Ij f 15. 191 ) induces an analogous 
splitting of the magnetic energy ( 15.101 ). i.e. 

Mim)^2 = E^ + Es, (5.22) 
with the orbital part Ej^ being found as 

El = -/II ■ i?ex. = ^eMex ■ j d¥ (r x kd{f)) , (5.23) 
and similarly for the spin part Eq, 

Es = -/Is ■ i?ex = ^e,#ex j (f X kp{f)^ . (5.24) 

Now it appears as a matter of course that when the RST correction Mi'^^c^ ( 15.221 ) is to 
be identified with the external magnetic energy Ei^t ( 15. (il ). the orbital part Ei^ ( 15.231 ) must 
be identified with the first part of E'int, i.e. 

and similarly the spin part Es f 15. 241 ) with the second part of Ei^t'- 

Es 2^#ex ■ j dV- (ij{f) ■ I ■ V'(r)) . (5.26) 
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Indeed, this claim f l5.25l ) - f 15. 261 ) can easily be verified: first consider the orbital part Ej^ 
f 15. 231 ) and observe that the drift current kd{r) reads in three - vector notation ( to be deduced 
from its four - vector version f 15. 151 ) ) 



(5.27) 



2Mc 

Substituting this into the definition of the orbital magnetic moment /II f 15. 201 ) and observing 
the Hermiticity of the angular momentum operator L yields immediately 

/7l = ^ [ d^f (i> ■ L ■ i/j) . (5.28) 



2Mc _ 

But with this result, the orbital magnetic energy Ei^ ( 15.231 ) is easily seen to coincide with 
its expected form f 15. 251 ) and thus the orbital part of the magnetic interaction energy E^^t 
( 15. (jl ) is in perfect agreement with RST. 

By a similar reasoning it is also possible to verify the expected correspondence of the 
spin part f 15. 261 ). albeit by means of a somewhat more subtle argument. For this purpose, 
observe first that the space part kp of the spin polarization current ^^^k^ f 15. 161 ) reduces to a 
three - curl if the stationary field configurations ( I3.1al ) - ( I3.1bl ) are considered: 

kp{r} = V X Psir) . (5.29) 

Here the spin polarization density Ps{r) is defined by 

Ps{r) = ^^^{r)-^-^{r) , (5.30) 

where the spin - operator 5* is introduced through 



l=^a (5.31) 



Now, when that polarization current kp{f) f l5.29l ) is inserted into the definition of the spin 
magnetic moment /is ( 15.211 ) , one ends up with the following final form ( by means of inte- 
grating by parts ): 



/?s = J d'f^ -S-^. (5.32) 

But clearly, with this result the RST spin magnetic energy E^ f 15. 241 ) is again identified 
with the spin magnetic part of the interaction energy i^int ( 15.61 ) and this verifies the claimed 
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correspondence ( 15.260 . 

Summarizing, the RST proposal f l5.8l ) for the external magnetic interaction energy is 
found to be in full agreement with the non - relativistic description of the (experimentally 
well-established) Zeeman effects; and this in turn may be taken as confirmation that also 
the internal magnetic interparticle interactions are correctly taken into account by the RST 
energy functional f 13. 391 ). 

VI. GROUND - STATE INTERACTION ENERGY 

For a demonstration of the results obtained so far, one may choose the simplest member 
of the para -system, which is the two -electron ground -state in the Coulomb field ( 13.91 ). 
Clearly, the arguments do apply also to all those excited states which own the same symme- 
try as the two particle ground- state, i.e. the highest possible symmetry. The simplification 
originates here from the fact that for the electrostatic approximation the mass eigenvalues 
Ma and the spatial parts of the wave functions ^"■'>R±(r) become identical, see equations 
( I3.15al ) - ( 13.181 ): and furthermore the time components Bo{r) and /io(0 of the exchange 
vector potential and exchange current /i^ do vanish, which is consistent with the Pois- 
son equation ( I3.7bl ). Furthermore, the identity of both mass eigenvalues implies that the 
exchange length om f 14. 201 ) becomes infinite which then annihilates the "electric" exchange 
field X (r ) f l4.21al ) ( -w X (r ) = ) . As a consequence there are no "electric" exchange cor- 
rections {m!i^^ = 0, see equation ( 14.241 ) ). and thus the "magnetic" exchange corrections 
(I1I2S1) simphfy to 

zi ■ Mi'^c" = Z2 ■ M^'^c^ = -P— [ dVf *(r) • Y{r) . (6.1) 

47ras J 

However, since the "magnetic" exchange field y(r) and its vector potential B{r) f l4.21bl ) 
do persist, they give rise to the emergence of exchange corrections for the para- system which 
are beyond the electrostatic approximation. The exchange vector potential B{r) obeys now 
the ordinary Poisson equation ( cf. equ. ( 14.331 ) ) 

AB{r) = -4:7iash*{r) . (6.2) 

with the corresponding simplification of the solution ( 14.361 ) 

B{f) =as f d'f'X^ ■ (6-3) 
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Observe also that the exchange three - current h{r) must become sourceless on behalf of the 



vanishing of its time component ho{r) ( 14.341 ) 

V-K*(f)=0. (6.4) 

Thus the Maxwellian exchange system ( I4.37an - f l4.37dn becomes truncated to its "mag- 
netic" part 

V ■ f(f) = (6.5a) 

V xY{r) = 47ras-h*{r) . (6.5b) 

Summarizing, the lowest - order corrections beyond the electrostatic approximation are 
described by the magnetostatic fields Ha{r) ( 14.181 ) and the "magnetic" exchange field Y{f) 
( 16.11 ) so that the total correction Ae!^^^ ( 14.521 ) becomes for the ground - state of the para- 
system 

a4-^) ^ -V^-Mi^^^^c^ (6.6) 



2^ 

a=l 



2 

Since this energy correction consists of the proper magnetostatic contributions (iVfi"^^ )and 
the "magnetic" exchange contributions ( Ma^^ ) , both subsystems have to be inspected now 
separately. 

A. Exchange Corrections Adjf^ 

Obviously, it is merely a technical question whether one prefers to compute the exchange 
masses Mi^^ in their original form f l4.9cl ) - f l4.9dn . being based upon the simultaneous use of 
the exchange current h(f) and the exchange vector potential B{r), or in the form f l4.42al ) - 
f l4.42bn which makes use of the currents alone ( a = 1,2) 



Za ■ M^^'c' = e JJ drd^r' ' , (6.7) 

or whether one prefers to deal with the above form ( 16.11 ) relying exclusively upon the "mag- 
netic" exchange field strenghts Y{r). In any case one must know explicitly the functional 
form of the exchange current h{r). Therefore one first inserts the general stationary form of 
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the wave functions ipai'f^, t) f l3.1al ) - f lS.lbH into the general definition of the exchange current 
h^{f) ( I3.8bl ) and then finds the following form for the three - current h(r): 



h{r) = — M+(r) ■ W;i^,ifi) . 



(6.8) 



Here the three -vector Wp depends on the spherical polar coordinates i?, (/? in the following 

way 

Wj^d, = - cos d (4 + iCy) + e'^P sin d , (6.9) 

where e^, e^, are the basis vectors due to a Cartesian parametrization ( x, 2; ) of Euclidian 
three -space. Futhermore, the radial function R_(_ has been defined in terms of the ansatz 
functions '^°''>R±{r) for the Pauli spinors ^'^'>(j)± as shown by equation f 13. 311 ). Since we are 
satisfied for the moment with the lowest -order approximation for the exchange masses 
we can resort to the non- relativistic approximation for the function ]R_|_ f 13. 321 ). 
With the exchange current h{r) ( 16.81 ) being at hand now, one can in the next step look for 
the solution B{r) f l6.3l ) of the Poisson equation ( 16.21 ). Clearly, the desired vector potential 
B{r) will have the same symmetry as its source h*{r), i.e. one tries the product ansatz 



B{r) =i-rB{r)-W.p{§,^) 



(6.10) 



and then one deduces the following differential equation for the radial ansatz function B[r) 
from the Poisson equation f l6.2l ): 

d^Bir) AdBir) Mi 



dr^ r dr 
The solution of this equation is easily worked out as 



as ■ 



(6.11) 



B(r) 



y 



Ob 



Ob 



+ 



-1 



Ob 



Ob 



■ exp ( —2 



Ob 



exp ( —2 



Ob 



- 1 



(6.12) 



Obviously, this solution decays as r at infinity ( r — 00 ) but is regular at the origin 
(r = 0) 



Bir) 



Ob 



+ 2- 

3 a-Q 



+ 



(6.13) 



The constant B^, is related to the normalization constant N^, f 13. 271 ) through (non-relativistic 
limit) 

B. = ^ (j^) ■ ^-"s i'^MNl) = ^alaj, . (6.14) 
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But once the exchange vector field B{r) is known, one can consider its curl Y(r) f l4.21bl ). 
which appears in the following form 

Yif) = -B{r) % + ley] + ■ f , (6.15) 

with the radial function B{r) being given by 

ABir) 

B(r)=r^^^ + 2B{r) . (6.16) 
dr 

When this result is used in order to compute the exchange corrections M^^^ ( 16.11 ). one is 
left after the the angular integration with the following radial problem (a = 1, 2 ): 

oo 

Ahr r r ~ 1 

Za-Mlf^^ = drr^ B{rf + 2B{rf . (6.17) 

3 as i L -I 



However, observing here the specific functional form of B[r) ( 16.161 ) and repeatedly inte- 
grating by parts yields 

oo 

Za ■ M[3^c^ j ^^B{r) ■ R+(r) , (6.18) 



where the differential equation for B{r) ( 16.111 ) has also been used. Now it is just this latter 
form f 16. 181 ) for the exhange corrections iVfi^'' which can also be recovered by starting from 
their original form f l4.9cn - f l4.9dn and using hereby the previously found results for the 
vector potential B{r) f 16.101 ) and the exchange current h{r) f l6.8l ). 

Clearly, the general equivalence of both forms ( I4.9cl ) - ( I4.9dl ) and ( 14.251 ) for the exchange 
corrections iVfi^^ has thus been exemplified merely in the lowest order ( beyond the electro- 
static approximation). But this can just be taken as a successfull consistency test of the 
applied approximation technique, which consists in a combination of linearizing the mag- 
netic interactions and additionally taking the non - relativistic limit. Indeed, the remaining 
radial integration of equation f 16. 181 ) can easily be done by use of the known functional forms 
of B{r) (Eni) and R+(r) (ISIS21) which yields 

z, ■ m['^c' = z, . M!f^c' = ^ (ze.«s)' .^ = ± (z^^asY E^r^ ■ (6-19) 

6 a-Q 15 

This result verifies the expectation that the "magnetic" interactions are typically smaller 
than the "electric" ones f 13. 501 ) by a factor [z^y^a^)'^ and therefore may be approximated 
here by their non -relativistic limit. But observe on the other hand that the magnetic 
corrections f 16. 191 ) vary as z^^ and therefore will become more important for the heavy 
atoms in comparison to their electric counterparts E)^ ( 13.501 ) raising only linearly with z^.^. 
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B. Relativistic Normalization 



It must be remarked that the value of the renormahzed "charges" Za f 13. 351 ) do not 
enter exphcitly the result f 16. 191 ) so that the values of the corresponding renormalization 
integrals f 13. 351 ) are not explicitly needed. But clearly, the renormahzed charges Za do enter 
the results in an implicit manner, namely through the proper relativistic normalization 



conditions |2l|. However, it is easy to demonstrate that for the presumed linearization 
f l4.10al ) - f l4.10bl ) the renormalization of the charges is trivial, i.e. one can identify: Za = 1. 
The reason for this is that the deviation of Za from unity is induced by the integral of the 



entanglement vector 



11 



20 



G, = ^^ [B-G:^ - B*^G,,] (6.20) 
over the time slices t = const which gives for the stationary field configurations 

z,-Z2 = - {Z2 ~Z2)= [ G^dS^ =^ [ d¥ Go(f) (6.21) 



<=const 

d^rB*(f) ■ Bir) 



27rasaM 

However, for the ground state the exchange length om becomes infinite ( om oo ) and 
thus the difference ( 16.211 ) between Za and Za vanishes: Za = Za = 1 , {a = 1,2). 

This coincidence of the charges Za and Za is not only a byproduct of the applied lineariza- 
tion, but for the ground state it has a deeper origin. Indeed, it is closely related to the 
symmetry of the two -particle ground- state. Observe here again that, when the normaliza- 
tion integral is done over a time slice ( t = const ) of the Minkowskian space - time, it is only 
the time component Go of the entanglement vector G^ ( 16.201 ) which is relevant: 

Go(f) = \B*{f) ■ X{r) - B{r) ■ X*{f)\ . (6.22) 

47ras L J 

This component however vanishes, because the "electric" exchange vector X(r) must be 

put to zero when its source ho{r) vanishes, see equation f l4.37an - f l4.37bl ). The latter fact, 

however, is a consequence of the symmetry of the ground - state ansatz f l3.1al ) - f l3.3f1 ). 

C. Magnetic Corrections Afi"*^ 

In order to complete the magnetic part AE^j^^'' ( 16.61 ) of the interaction energy, we have 
to consider the magnetostatic mass corrections AijT'^ ( I4.7al ) - ( I4.7bl ). This will be done by 
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explicit computation of the magnetostatic vector potentials Aa{r) f l4.29bn . however through 
solving directly the associated Poisson equation f l4.27bn . instead of computing the integral 
f l4.29bn . The reason for this is that both magnetostatic potentials Aa{f) (a = 1,2) are of 
a similar functional form as the "magnetic" exchange potential B{f) f 16. 101 ). so that one 
can take over the functional form of the corresponding solution for the desired potentials 
Aa{r). Clearly, this similarity of the vector potentials B(r) and Aa{r) has a deeper geometric 
meaning to be discussed readily. 

The close relationship between all three potentials B{r), Aa{r) {a = 1,2) is recognized 
most immediatly by demonstrating that the Poisson equations for the magnetic potentials 
Aa{r) f l4.27bn are efiictively the same as for the exchange potential B{r) ( 16.20 . i.e. the 
right - hand sides of both Poisson equations effectively do agree. This does not mean that the 
currents ka{r) are identical to the exchange current h{f) ( 16.81 ). but when one computes the 
three- currents ka{r) ( KlSal ) by means of the stationary field configurations ( Kllal ) - ( I3.1bl ). 
one finds the following form: 

kaif^ = kair)-Vpi^,ip) . (6.23) 

Here the vector field Vp is given by 

Vp{-^, if) = sin-i? ■ [— sinyj ■ + cosy? ■ Cy] = sin?? ■ (6.24) 

and thus is to be conceived as the magnetic analogue of its exchange counterpart Wpi^^ ip) 
( 16.91 ). Furthermore, the scalar prefactors ka{r) in equation ( 16.231 ) are found to be of the 
following form 

ki{r) = — Wi?+(r) ■ Wi?_(r) (6.25a) 

271 

k2{r) = - — (2)i?+(r) ■ (2)i?_(r) . (6.25b) 
27r 

However, since the radial functions of both particles do coincide for the ground- state (i.e. 
(i)i?+(r) = (2)i?+(r) , (i)i?_(r) = (2)i?_(r) ) one has 

ki{r) = -k2{r) = -^M+(r) , (6.26) 
An 

and therefore both single- particle currents ka{'r) ( 16.231 ) flow in opposite directions around 
the 2; -axis (see fig.l) 

fci(r) = -hir) = k{f) . (6.27) 
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FIG. 1: Isotropy of the Ground - State 

Each coordinate axis is encircled by the same kind of circular flow k f l6.27l ). ^ f l6.28bl ') and ff 
f l6.28cl ) which yields the equipartition of the magnetic interaction energy Ae!^^^ (EM) into the 
identical contributions AE^^'^ (EM), A-E^f ^ ( \(\Aia\ ) and Ae!^'' f l6.41bn . 

Now recall the expectation that the ground - state symmetry has to be the highest possible 
one, so that the 2;- axis cannot be singled out in any way as compared to the x- and ?/- axis. 
Indeed, this is really the case if one takes also into account the exchange current h{f) f l6.8l ) 
which splits up into its real and imaginary parts C,{f) and f/(r) as follows: 

h(f^ = ^(r) + if]{r) (6.28a) 
-• M_|_(r) 

^(r) = — [— cos ■ e„ + sin sin ip ■ Cz] (6.28b) 

47r 

vi^ — — [~ '^os ■ Cx + sin^ cos (f ■ Cz] . (6.28c) 

4:71 

Evidently, this result says now that all three axes are encircled by the same type of flow, 
namely the x-axis by ^(r) f l6.28bn . the y-axis by ff{f) f l6.28cl ) and the z-axis by k{r) 
( 16.271 ): see fig.l. Here it must be stressed that this isotropic geometry for the ground -state 
could only be attained by completing the magnetostatic potentials Aa{'r) by the exchange 
potential B{f) into a triplet of vector potentials. Since B{f) owes its existence to the use 
of the non-abelian group U{2) instead of the abelian U{1), which is generally used in the 
conventional approach to the electromagnetic interactions ( classical and quantum ) , the 



56 



"magnetic" exchange interactions are now seen to equip the RST with a truly non- abelian 
character! 

For the exphcit construction of the magnetostatic potentials Aa{f), one tries an ansatz 
which has the same symmetry as their sources ka{r) f 16.231 ). quite similarly as it was done 
for the exchange potential B{f) f l6.iun . i.e. one puts 

Ai{r) = rAi(r) ■ Vpi^, if) (6.29a) 
Mr) = rA^ir) ■ Vp{^, ^) . (6.29b) 

Inserting this into the Poisson equations f l4.27bn yields the following differential equations 
for the radial ansatz functions Aa{r): 

j!diM ^ i jMl ^ 2„."'^^''-' ■ "'^-<'-) . 4.t,(.) (6.30a) 

ar^ r ar r 

+ l^fl^ = 2as ' = Aixk^ir) . (6.30b) 

dr^ r ar r 

Since for the ground - state symmetry the radial parts of both wave functions '4)a{r) do agree 
(i.e. (^)i?+(r) = (2)i?^(r) , (i)i?_(r) = (%„(r) ), one arrives at 

dMi(r) , 4 dA,(r) M+(r) 

— ^ = «s (6.31a) 

dr^ r dr r 

dM2(r) , 4 dA2(r) M+(r) 

— \ 1 = -as , (6.31b) 

dr^ r dr r 

which is effectively the same equation as for the radial exchange function B{r) f 16.110 . 
Therefore, one can directly take over the solutions as 

Ai(r) = -A2{r) = B{r) , (6.32) 

where the radial function B{r) has already been specified by equation ( 16.121 ). 

Clearly, the identification ( 16.321 ) of all radial functions Aa{r) and B{r) is a further in- 
dication of the ground -state isotropy and it is instructive to elucidate this effect also from 
another point of view. For this purpose, reconsider the source equations f l2.40an - f l2.40dl ) 
and insert therein the general shapes of the currents ka{r) f 16. 231 ) and h{f) ( 16.81 ). together 
with the vector potentials Aa{r) f l6.29an - f l6.29bl ) and B{r) f 16. 101 ). The right-hand sides 
of the first two equations ( I2.40al ) - ( I2.40bl ) then turn out to be zero 
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(6.33) 



and therefore the currents ka{r) are found to be sourceless 

V ■ = . (6.34) 

Observe that this property is also shared by the exchange current h{r), cf. f l6.4l ). The 
remaining two equations ( I2.40cl ) - ( I2.40dl ) yield the relation 

Ai(r)-A2(r) M+(r) 
^i^) = 1 I \ — n~\ 'a — • (6.35) 

Now recall here that, because of the ground- state symmetry, the denominator becomes on 
account of equation ( 16.261 ) 

ki{r)-k2{r) = ^^ (6.36) 
Zn 

and furthermore the numerator appears by virtue of equation f 16. 321 ) as 

Ai{r) - A2ir) =2B{r) (6.37) 

so that the relation f 16. 351 ) degenerates to an identity. Therefore the source equations 
f l2.40an - f l2.40dn are automatically satisfied when all RST objects | A;a(r), h{r), Aa{r), B{r) | 
own the symmetry properties of the ground - state, irrespective of the precise functional form 
of the radial parts '•"^i?±(r) of the wave functions! 

With these prerequisites at hand, we are now able to test a further consequence of the 
ground -state isotropy: if it is true that, on account of the isotropic geometry, all three 
spatial directions are equivalent, then the corresponding currents k{f) f 16. 2 71 ). ^(r) f l6.28bl ) 
and r/(r) ( I6.28cl ) have to contribute equal parts to the magnetic energy Ai?,^™'^'' f l6.6l ). see 
fig. 1: 

a4™'^ = a4"^ + 4'^ + 4'^ • (6-38) 
Now the contribution of the 2; -direction is given by the currents ka{r), cf. ( 14.391 ) 



a4^) ^ ^ (-^1 ■ ^t^c^ + ^2 ■ Mt^^) =e^ jj rfV-rfV--'^^^-^ . (6.39) 

Next, the contributions of the a; -and directions are given by the "magnetic" exchange 
energy Za ■ Mjf^c^, see equations ( I4.42al ) - ( I4.42bl ) for infinite exchange length {au ^ oo): 

a4") + A^) ^ ^ (^1 ■ + Z2 ■ M^'^c'^ =e^ II d^fd'f' , (6.40) 
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i.e. for either direction separately 

= e^ jj d'fd^f' ^^^'J^^^l^ (6.41a) 
A^) =e' JJ d'fd:'r' "^^^'J^"^^^ ■ (6.41b) 

But since the three currents k{r),^{f),f]{r) differ only by their spatial orientation, not by 
their intrinsic pattern of flux lines and strength, all three contributions are identical, cf. 

(Eini) 

Ae!^^^ = AE^ + A^) + A^) = 1 {z^^asf " ^ = ^ (^e.«s)' ■ • (6.42) 

This is the definitive result for the magnetic corrections; and can now to be used in 
order to check our preliminary estimate of the magnitude of these magnetic contributions 
which has been made within the framework of the electrostatic approximation ( see end of 
sect. III). Namely for small values of the coupling constant ( Zex^s 1); the preliminary 
picture of the magnetic interactions did identify the electrostatic RST prediction Aii^p^grp 
with the non - relativistic interaction energy E)^ , see equation f 13. 491 ). On the other hand 
the discrepancy between these RST predictions Aii^^grp and the experimental data Aii^exp 
was attributed to the magnetic energy A^™^"*, see equation f 13. 5 71 ). Therefore the limit 
value /o of the geometric factor f 13. 581 ) for [z^^asf <^ 1 appears now as 

/o = 7 ^ • ^f- = T (6.43) 



which appears to be in acceptable agreement with the values of as presented in table I. 
Clearly, for a more precise test of the limit values ^3.521 ) and /q f 16. 431 ) one would have 
to extend table I to smaller values of the coupling parameter ^ex«S) possibly up to neutral 
helium ( separate paper ). 



D. Theory vs. Experiment 



Apart from such an internal consistency test for the RST predictions, it is of course highly 
^cresting to compare 
magnetic contributions 



interesting to compare the complete RST predictions Ai?^^"*, i.e the sum of electric and 



A^s^) = AE^T + a4-^) , (6.44) 
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Element 

(^ex) 


[12] 


(Km) 




Relativistic MBPT 


Rel. all order 
MBPT [19] 


MCDF 

M-m 


Unified 
[13] 


^Ecxp 


Ge ( 32 ) 


562,5±1,6 


564,9 


-0,42% 


561,9 


562,1 


562,1 


562,1 


Xe ( 54 ) 


1027,2±3,5 


1031 


-0,37% 


1028,1 


1028,4 


1028,2 


1028,8 


Dy (66) 


1341,6±4,3 


1336 


0,41% 


1336,6 


1337,2 


1336,5 


1338,2 


W (74) 


1568±15 


1570 


-0,11% 


1574,6 


1574,8 


1574,6 


1576,6 


Bi ( 83 ) 


1876±14 


1868 


0,43% 


1882,7 




1880,8 


1886,3 



TABLE II: Comparison of the theoretical predictions for the ground -state interaction energy 
with the experimental data All energies are measured in [eV]. Already after the inclusion of 
the magnetic interactions in the lowest - order approximation, the RST predictions AE"^^^ ( third 
column ) meet with the corresponding predictions of the other theoretical approaches ( right half 
of the table ). 



to both the experimental data and to other theoretical predictions, such as the relativistic 
1/Z expansion the multiconfiguration Dirac-Fock method (MCDF) ^|- ^| and rel- 
ativistic many -bo dy p ertubation theory (MBPT) all -order technique for 
relativistic MBPT [19]. The predictions of these four theoretical approaches are collected 



in table II together with the experimental data [12] in order to oppose them to the present 
RST predictions AE^^^t (EM)- As it can easily be seen from a comparison of table I and 
table II, the RST predictions appear now to be very satisfying: after the magnetic interac- 
tions have been included into the RST results to yield the interaction energy AE^^^ ( 16.441 ). 
third column of table II, the deviation from the experimental data Ai^cxp ( second column 
of table 2 ) has decreased to less than a half percent which does no longer depend upon the 
nuclear charge Zex- This means that the discrepancy between the experimental values A£^exp 
and the electrostatic RST predictions Aii^^grp, ranging from 1,7% for germanium( 32 ) up to 
11,5% for bismuth ( 83 ) ( see table I ), is actually caused by the magnetic interactions with a 
remaining uncertainty of roughly 0,4%. This supports the RST picture of the simultaneous 
action of electric and magnetic forces in the electronic orbits of an atom. 

A further, very satisfying feature of the RST results AE^^^^ ( 16.441 ) refers to the com- 
parison to the other four theoretical approaches, see the last four columns of table II. Here 
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the RST predictions appear to be of the same order of precision as the 



corresponding pre- 



12| . Summarizing, 



dictions of the other four theoretical approaches, see table III of ref. 
one can therefore judge that the RST predictions can compete succesfuUy with these other 
approaches, even when the magnetic interactions are treated in the lowest-order approxi- 
mation! Clearly, in the next step one will study the higher -order approximations of the 
magnetic RST interactions, which may be expected to lead even closer to the experimental 
data. 
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